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ABSTRACT 

We present a spectro-photometric survey of 2522 extragalactic globular clusters (GCs) around 
twelve early-type galaxies, nine of which have not been published previously. Combining 
space-based and multi-colour wide field ground-based imaging, with spectra from the Keck 
DEIMOS instrument, we obtain an average of 160 GC radial velocities per galaxy, with a 
high velocity precision of ~ 15 km s _1 per GC. After studying the photometric properties 
of the GC systems, such as their spatial and colour distributions, we focus on the kinemat- 
ics of metal-poor (blue) and metal-rich (red) GC subpopulations to an average distance of 
~ 8 effective radii from the galaxy centre. 

Our results show that for some systems the bimodality in GC colour is also present in 
GC kinematics. The kinematics of the red GC subpopulations are strongly coupled with the 
host galaxy stellar kinematics. The blue GC subpopulations are more dominated by random 
motions, especially in the outer regions, and decoupled from the red GCs. Peculiar GC kine- 
matic profiles are seen in some galaxies: the blue GCs in NGC 821 rotate along the galaxy 
minor axis, whereas the GC system of the lenticular galaxy NGC 7457 appears to be strongly 
rotation supported in the outer region. 

We supplement our galaxy sample with data from the literature and carry out a number 
of tests to study the kinematic differences between the two GC subpopulations. We confirm 
that the GC kinematics are coupled with the host galaxy properties and find that the velocity 
kurtosis and the slope of their velocity dispersion profiles is different between the two GC 
subpopulations in more massive galaxies. 

Key words: galaxies:star clusters - galaxies:evolution - galaxies: kinematics and dynamics 
-GCs 



1 INTRODUCTION 

The reconstruction of the evolutionary history of galaxies in the lo- 
cal universe requires a comprehensive knowledge of their chemo- 
dynamic properties at all scale s. Along these lines, the ATLAS 3D 
project dCappellarietal.ll201ll) has carried out a rich survey of 
galaxies employing an integral field unit (IFU) technique to map 
out the chemo-dynamics of galaxies within 1 effective radius 
(Res)- However, such a radius encloses < 10% of the total galaxy 
mass (baryonic + dark) and hence it may not be representative of 



the overall galaxy (e.g., IProctor et al]|2009l) . In contrast, beyond 
this radius (7? c ff > 1), the halo of the galaxy contains precious dy- 
namical signatures of early merging or earl y collapse events that 
migh t have eventually built up the galaxy dHopkins etai]|2009[ 
I2O10I : lOser et al.|[20lol ; iHoffman et afll20ld) . Some recent exam- 
ples of the dynamical studies of the outskirts of galaxies beyond 
the Local Grou p has mainly involv ed a few deep long-slit ob- 
servations (e . g.,ICoccato et alj|20ld). the pi oneering technique of 
IProctor et all ( 120091) and lNorris et al.1 d2008h that probed the two- 
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dimensional stellar kinematics to 3_B e ff , and extragalactic plane- 
tary nebulae (PNe, ICoccato et ai]|2009I) that can probe the galaxy 
potential up to ~ 10i? o ft . 

Often labeled the "fossil record" of galaxy formation, GCs 
have been used to study galaxy haloes. GC systems have been 
found in galaxies of all morphological type and they extend 
beyond the detectable galaxy diffuse light, probing galactocen- 
tric distances > 10R e B where the are only few other gravita- 
tional constrains dRhode & Zepfll200ll . 12004 ; iDirsch et all 120031 ; 



2000; Forbes et al 



iTamura et al.l2006 ; Forbes et alj201 ll) Moreover, their old ages (> 



10 Gyr.lKissler-Patig et al. | ll998l;ICohen et al.lll998l;lBeaslev et al 



20021 ; iBeaslev et al 



2001; Schroder et al 



2002T 



Brodie & Larseri 



|2004| ; IStrader et al.l l2005h suggest they have 



survived violent merging events, preserving the chemo-dynamical 
record of their parent galaxies. 

Moreover, the well-studied dichot omy observed in the 
colour distribution of most GC sy stem s 1 Zepf & Ashman! 19931: 



colour distribution ot most OL, systems (Zepl& Ashman 1993; 
Ostrov et al.ll993l;IWhitmore et al.ll99l;lGebhardt & Kissler-Patig 



1999 1: lLarsen et ail l200ll; iKundu & Whitmord 1200 ll ; IPeng et al 



2006 ; Sinnott et al .1 201fJ) is thought to stem from different 
formation mechanisms that shaped the underlying host g alaxy 
jAshman&ZepJll992l ;l Forbes et al. Hl997l ; IC6teetai]|l998l) . Al- 
though the reality of the GC metallicity bimodality has been re- 
cently called into questio n as the result of a strongly non-linear 
color -metallicity relation dYoon et al.ll2006l I20ID ; iBlakeslee et all 
l20ld) . observations have shown th at physical dissimilarities ex- 
ist among these two subpopulations 1 Cote 1999; Brodie & Strader 
20061; IPeng et al.l 120061 ; IChies-Santos et al] l201ll ; iForbes et al 



20111) . For instance, the t wo GC subpopulations are f ound to have 



different physical sizes dKundu & Whitmord l200ll ; | Jorda n et al .1 
120051 ; [Masters et al1l2010t) and diverse spatial distributions around 
the host galaxy, with the metal-rich (red) GCs more centrally 
concentrated than the metal-poor (blue) GCs dGeisler et al.lll99q ; 



Ashman & Zepf 1998; Brodie & Strader 2006; Bassino et al. 2006; 



Faifer et al.l201ll ; IStrader et al.ll201ll ; lForbes et alj|2012l) . 



Current spectroscopic studies of GC systems have shown that 
the kinematics (e.g. rotation directions and rotation amplitudes) of 
the two GC subpopulations are somewhat diverse. The kinemat- 
ics of the red GCs is usually akin to that of the host galaxy stars 
dSchuberth et alj20ld;ls"trader et alj201 lh . perhaps due to a similar 
formation history 1 Shapiro et alj|20ip|) . Whereas, the velocity dis- 
persi on of the blue G Cs is typically larger than that of the red GCs 
(e.g.. iLee et alj|2008h . Also intriguing is the fact that rotation has 
been detected for both the blue and the red GCs, regardless the mass 
or mo rphology of the host galaxy dFoster et alj|201ll ; lArnold et al .1 
l20Tlh . 

Interpreting this variety of GC kinematics in the context of 
galaxy and GC formation has been limited due to the low num- 
ber of galaxies with large GC radial velocity datasets. To date, this 
set includes only a dozen G C systems, most of which are nearby 
very massive ellipticals (see ILee et alj|201ol for a summary). On 
the other side, numerical simulations on this front have mainly fo- 
cused on the origin of GC metallicity and colour bimodality (e.g. , 



Weil & Pudritzll200ll ; iKravtsov & Gnedinll2005l ; lYoon et alj|2006l 
20111) . rather than on GC kinemat ic properties (e.g., iBekki et al .1 



2005ll2008l ; |Prieto & Gnedinll2008l) . 



This scenario has left open several questions regarding the 
kinematics of GC systems: does the colour bimodality also im- 
ply kinematic bimodality? If so, do the kinematical differences be- 
tween the blue and red GC subpopulations found in the most mas- 
sive ellipticals also hold for ~ L* galaxies over the whole early- 
type sequence of the Hubble diagram? Furthermore, do the blue 



and red GC subpopulations rotate faster in the outer reg ions, as pre- 
dicted in a formation in a disk-disk merging scenario l lBekki et"al] 
120051) ? Can GC kinematics contribut e to our understa nding of the 
formation of lenticular galaxies (e.g jBarr et alj|2007h ? 

We have been carrying out a project named SLUGGsO] to in- 
vestigate the GC systems in external galaxies (Brodie et al. 2012, 
in preparation). SLUGGS is the SAGES Legacy Unifying Globu- 
lars and Galaxies Survey, where SAGES is the Study of the As- 
trophysics of Globular Clusters in Extragalactic Systems. This sur- 
vey exploits the combination of Subaru/Suprime-Cam wide-field 
imaging with spectra from the Keck/DEIMOS multi-object spec- 
trograph. The results released so far have shown that the wide- 
field i maging can give clues about assemb ly history of the host 
galaxy dBlom et alfeoi ll ; lForbes et al.l20l"ll) . If combined with the 
high velocity resolution of DEIMOS, this dataset can unravel, at 
the same time, the dynamics and the metallicity of t he field stars 
dFoster et al.l2009l:|Proctor et al.l2009l) and of the GCs dFoster et all 
1201 ll ; lArnold et alj|201 lh deep into the galax y halo, as well as a 



I2U1 11 : Arnold et al. 2U1 JJ) deep into tne galax y nalo, as well as a 
giving compelling view of gal axy dynamics dStrader et alj|201 H : 
iRomanowskv et alj2009l.l201 lh . 

In this work, we aim to study the GC kinematics for an un- 
precedented large sample of early-type galaxies (ETGs) with high- 
quality data. We investigate the global kinematics of blue and red 
subpopulations to study how they, and the underlying galaxy itself, 
formed. We also supplement our galaxy sample with literature data, 
and we compare the properties of this large sample with the existing 
numerical predictions. 

The plan of this paper is as follows. From Section 1 to Section 
5 we describe the reduction and the analysis of both the photometric 
and the spectroscopic data. In Section[5] we briefly discuss the sig- 
nificant findings for each galaxy. In Section[6] we give an overview 
on the current state of the GC formation models in order to com- 
pare their predictions to our generic results discussed in Section]?] 
In Section[8]we supplement our galaxy sample with data from the 
literature. In Section|9]we analyse the enlarged GC system dataset 
(our data plus literature) and discuss the results in SectionQT)] The 
summary of the paper is given in SectionQT] 



2 THE SAMPLE 

In this paper we discuss a subset sample from our survey. This 
includes nine new galaxies, in addition to other three galaxies 
already published: NGC 4494 dFoster et alj 1201 ll). N GC 3115 
d Arnold et al.l20 111) and NGC 4486 dStrader et al J201 ll) . The anal- 
ysis and the specific results for these three galaxies have been ex- 
tensively discussed in the respective papers. Therefore, their overall 
results will be discussed together with the other nine starting from 
Section|7j 

The physical characteristics of the twelve galaxies are listed in 
Table Q] with their optical images shown in Figure Q] This galaxy 
sample extends the study of extragalactic GC systems into a new 
regime, because it is representative of a wide range of luminosity, 
morphological type (from lenticulars to giant ellipticals) and envi- 
ronment (from field to clusters), with a velocity resolution three 
times better than typical previous studies. This improvement is 
shown in Figure [2] in which we compare the intrinsic properties of 
our dataset with previous GC studies (that have employed various 
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Galaxy ID 


Hubble 


Vsys 


Z? c ff 


(to - M) 


D 


M K 


A K 


PA K 


(b/a) K 




Type 


[kms" 1 ] 


[arcsec] 


[mag] 


[Mpc] 


[mag] 


[mag] 


[degree] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


NGC 0821 


E6 


1718 


51 


31.85 


23.4 


-24.0 


0.040 


31 


0.62 


NGC 1400 


SAO 


558 


31 


31.05 


26.8 


-23.8 


0.024 


29 


0.90 


NGC 1407 


E0 


1779 


72 


32.30 


26.8 


-25.4 


0.025 


70 


0.95 


NGC 2768 


E6 


1353 


93 


31.69 


21.8 


-24.7 


0.016 


91 


0.46 


NGC 3377 


E6 


690 


46 


30.19 


10.9 


—22.7 


0.013 


46 


0.58 




E2 


620 


34 


30.97 


15.6 


—23.8 


0.010 


219 


0.91 




Co 




DU 




no q 


— zo.z 


U.UUo 




U. f 4 


NGC 5846 


E0 


1712 


61 


31.92 


24.2 


-25.0 


0.020 


233 


0.92 


NGC 7457 


SO 


844 


52 


30.55 


12.9 


-22.4 


0.019 


125 


0.54 


NGC 3115 


SO 


663 


85 


29.87 


9.4 


-24.0 


0.017 


43 


0.45 


NGC 4486 


E0 


1284 


81 


31.18 


17.2 


-25.3 


0.008 


151 


0.86 


NGC 4494 


El 


1344 


53 


31.10 


16.6 


-24.2 


0.008 


173 


0.87 



Table 1. General properties of our galaxy sample. T he galaxy name ( 1 ) and Hubble Type (2) are from the NED database. T he galaxy systemic ve locity (3) 
and ellipticity correcte d effective radius (4) are from Cappellari et alj i201 ll) . otherwise from NED and the RC3 catalogue (d e Vaucouleurs 1 99Tb if not in 
ICappellari et alj 1 201 1|). The d istance modulus (5) and the respective distance in Megaparsec (6) are from lTqnry et alj l200l|) with a —0.06 correction as 
advocated bv lMeiTt*afU2007l) , respectively. If the galaxy is in the ACS Virgo Surv ey, we use the dista nces from lMei*et*alli2*007l) . We assume that NGC 1407 
and NGC 1400 lie at the same distance, computed as the average of the respective iTonrv et all 1200 "fj ) distances. The K band absolute magnitude (7) is from 
2MASS apparent magnitude at the distances given in column 6 and corrected for the fore ground Galactic extinction given in column 8 (NED database). The 
photometric position angle (9) and axis ratio (10) are from 2MASS ISkrutskie et alj200rj) . The last three galaxies have been analysed in separate papers (see 
text). 




NGC 3377. ' 'J... NGC 4278' ' / • NGC 4365 : ' o ■ .NGC 5846 



♦ 




Figure 1. Digitalized Sky Survey (DSS) images of our galaxy sample. Red circles and black ellipses represent the spectroscopically confirmed GCs found in 
this work and the 2, 4, 6 i? c ff schematic isophotes corrected for the galaxy ellipticity respectively. The black line on the bottom-right spans 5 arcmin in length. 
North up and the East on the left. The elliptical galaxy south-west of NGC 1407 is NGC 1400 whose GC system is shown separately in this figure. The last 
three galaxies have been analysed in separate papers. 



© 2012 RAS, MNRAS 000,[[|-?? 



4 Pota et al. 



60 

40 



3- 20 



,2 2 



....... 

Size oc -^N~ 


ill 


....... 

. ■ : 




■ 

. ■ 


-■ : 


o 

o 






— 1 — 1 — 1 — 1 — 1 — 1 — 1 — 
Size oc 1 j 7 (SV) 

- o 


1 — 1 — | — 1 — h- 


H — 1 — | — 1 — 1 — 1 — 1 — 

o'o . - 

° 

o ■ 


o 

, , , , i , , 


, i , , 


■ 

I , 1 I , , I 


-23 


-24 
M K [mag] 


-25 -: 



Figure 2. Accuracy of GC radial velocity measurements. Open circles and 
open boxes represent our GC data and literature data, respectively. Liter- 
ature data will be discussed in Section [8] Top panel. The median velocity 
uncertainty (SV) per GC system are shown as a function of absolute mag- 
nitude in the K band. Symbol sizes are oc V~N, where N is the number 
of spectroscopically confirmed GCs per galaxy. Bottom panel. Log JV is 
shown as a function of the absolute magnitude in the K band, where sym- 
bol sizes are oc 1/ (SV). Our data extend the range of galaxy luminosity 
probed with three times better velocity accuracy. 



instruments including VLT/FLAMES, VLTVFORS2, Keck/LRIS or 
Gemini/GMOS). 



3 PHOTOMETRIC OBSERVATIONS AND DATA 
ANALYSIS 

3.1 Subaru data 

Multi-band photometric observations were carried out using 
the Subaru/Suprime-Cam instrument with a fie ld of view of 
34 x 27 arcmin 2 and a pixel scale of 0.202 arcsec l lMivazaki et al .1 
12002b . The galaxy sample was observed in the period between 
2005 and 2010 with a standard Sloan Digital Sky Survey (SDSS) 
gri filter set. NGC 2768 a nd NGC 4278 we re both downloaded 
from the SMOKA archive dBaba et al]|2002l) and were observed 
with a Rciz and BVI filter set, respectively. For NGC 1407 and 
NGC 43 65 we present exis ting Supri me-Cam photometry pub- 
lished in ISpitler et al.1 d2012l) (see also Romanowskv et alll200°j) 
and iBlom et al.l d2011 ), respectively. For these two galaxies both 
the imaging reduction and the catalogue extraction were performed 
with the methodology described in this paper. NGC 7457 is the 
only galaxy for which no Subaru imaging is available, and there- 
fo re we use WIYN/M inimosaic imaging in BVR filters presented 
in lHargis et al.l J201 ll) and we refer to this paper for a description 
of the data reduction. 

In Table [2] the imaging observations are summarized. The 
overall seeing conditions were mainly sub-arcsec. For NGC 821 
(g band), NGC 5846 (g band) and NGC 4278 (B band) the data 
suffer from cloudy conditions and poor (> 1 arcsec) seeing. 

Subaru raw images were processed using the SDFRED data 



pipeline dOuchi et alj|2004j) that yields standard flat field corrected 
images for each of the three filters. Photometric point source cat- 
alogues were extracted using standard IRAF/Daophot aperture 
photometry routines. We summarise h ere the main steps o f the data 
reduction and we refer to Section 3 of IBlom et al .1 d201lh for a de- 
tailed description of the method. 

We obtain a raw list of object positions by running 
IRAF/Daof ind on galaxy subtracted images in order to opti- 
mise the finding algorithm. The extraction threshold was typically 
set between 2 and 4 times the background depending on the filter 
and on the seeing conditions. The galaxy light was modelled with 
I RAF /Ell ipse set to allow the position angle and ellipticity to 
vary. Next, we perform aperture photometry using I RAF /Phot 
on the preselected objects for a certain number of circular aper- 
tures from 1 up to 15 pixels (equivalent to ~ 0.2 arcsec to 3 arc- 
sec for the Suprime-Cam pixel scale). The extraction radius was 
chosen in order to maximise the signal of the source and minimise 
the sky contribution. The extracted magnitude was corrected for the 
computed aperture correction using I RAF / Mkapfile. Photomet- 
ric zeropoints were estimated by boot-strapping the Suprime-Cam 
photometry to the Sloan Digital Sky Survey (SDSS) DR7 photo- 
metric system dAbazaiian et al.l200°j) using the brightest objects in 
common between the two datasets (typically with 17 < i < 21). 
If not in SDSS. the zeropoints were calibrated using the flux from 
standard stars observed over the same night. Finally, we use the red- 
dening given in Table Q] and the conversion table of ISchlegel et al.l 
dl998 ) to derive the Galactic extinction correction in our photo- 
metric bands. Hereafter, all magnitudes and colours are extinction 
corrected. 



3.2 HST data 

We use Hubble Space Telescope (HST) archive images from the 
Advanced Camera for Surveys (ACS) and Wide Field Planetary 
Camera 2 (WFPC2) to improve the quality of the photometric se- 
lection in the central regions of our galaxies. For most of them, we 
exploit existing photometric GC catalogues and we refer to the fol- 
lowing a uthors for a detailed description o f the data reduction and 
analysis: Ispitler et all d200^) for NGC 821 [Forbes et alJd2006t) for 
NGC 1407 and NGC 1400, iForbes et alj dl996t) for NGC 5 846, 
IChomiuk et all d2008l) for NGC 7457 and lBlom et ail d201lh for 
NGC 4365. 

We obtained, from the Hubble Legacy Archive, new HST/ ACS 
imaging for NGC 3377, NGC 2768 and NGC 4278, respectively. 
The ACS camera has a pixel scale of 0.05 arcsec and a field of view 
of 3.36 x 3.36 arcmin 2 . 

The NGC 3377 imaging consists of one pointing in F475W 
(~ Sloan g) and F850LP (~ Sloan z) filters and it was observed 
as part of the HST project ID 10554. NGC 2768 (ID 9353) was 
imaged in F435W, F555W, F814W filters, equivalent to a BVI 
configuration, respectively. Finally, the NGC 4278 (ID 10835) data 
consists of four pointings in F475W and F850LP filters that probe 
the galaxy up to ~ 6 arcmin from the centre (Usher et al. 2012, in 
preparation). 

The HST imaging was reduced and analysed using a custom 
built pipeline to find point-like sources and measure their magni- 
tudes and half light radii. For details on the methods used by the 
pi peline including poi nt sprea d function d e termin ation, we refer 
to IStrader etal] d2006t) and to ISpitler et all d2006l) . The extracted 
magnitudes and sizes for GCs in N GC 3377 and NGC 4278 w ere 
compared with those published by IChies-Santos et af] d201 lh for 
objects in the NGC 3377 pointing and in the two NGC 4278 point- 
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Figure 3. Colour-colour diagram of objects around NGC 1407. Grey and 
black points represent all the sources detected in the Suprime-Cam field-of 
view and all the photometrically selected GCs with i < 25 and i < 24, 
respectively. Spectroscopically confirmed sources are shown as red points 
(GCs), green boxes (Galactic stars) and orange triangles (galaxies) respec- 
tively. 



ings. Both magnitudes and sizes show good agreement without any 
evidence of statistically significant offset from the published data. 



3.3 Photometric GC selection 

Once the point-source catalogues have been constructed, they are 
photometrically selected to avoid contamination, such as unre- 
solved galaxies and Galactic stars. As far as our new Subaru data 
are concerned, this process involves the following steps: (i) selec- 
tion in a colour-colour space; (ii) selection in a colour-magnitude 
space; (iii) a visual check. The GC select ion in the ground-b ased 
imaging for NGC 7457 was carried out in lHargis et all d201ll) and 
we refer to this paper for a detailed description of their selection 
method. The steps listed above are subjected to variations in the 
procedures for the DEIMOS mask design. 

GCs are known to populate a specifi c area of colour- 



colour diagrams (e.g., 



Rhode & Zepl 1200 ll ; iFaifer et all 1201 ll ; 



IChies-Santos et all 1201 1 ). With a gri filter set, this is enclosed 
within 0.4 < (g - i) < 1.5, < (r - i) < 0.6, 0.3 < (g - r) < 0.9 
where these boundaries run diagonally to the colour axes, as shown 
in Figure[3]for the galaxy NGC 1407. To take into account the de- 
pendence of the colour boundaries on the quality of the data, we 
flag as GC candid ates all the objects deviating by less than 2<r from 
these boundaries JSpitler et alj|2008h . 

Next, we apply a cut on the i band magnitude (/ band 
for NGC 4278). Given the ongoing debate regarding the uncer- 
tain separation be t ween GCs and ul tra compact dwarfs (UCDs, 
iMieske et al.ll2006| ; lBrodie et alj201 l|), we decided to set the upper 
brightness magnitude at Kh ~ —11.6 (Mi ~ —12), one magni- 
tude brighter than the integrated magnitude of oj Cen, the brightest 
GC in the Milky Way. Nevertheless, in some cases we relax this 
criterion in order to include spectroscopically confirmed GCs that 
have magnitudes brighter than the set threshold. The separation be- 



Galaxy ID 


Obs date 


Filters 


Exp. time 


Seeing 








[sec] 


[arcsec] 


NGC 0821 


2008 Nov. 


gri 


960, 350, 220 


1.2,1.0,0.8 


NGC 1400 


2008 Nov. 


gri 


3240, 3600, 10800 


0.6,0.6,0.6 


NGC 1407 


2008 Nov. 


gri 


3240, 3600, 10800 


0.6,0.6,0.6 




iUuJ IVXcU. 




Ron 1 no 91 n 

UUU, -LUU, ilU 


u. u, u. u, u. u 


NGC 3377 


2008 Nov. 


gri 


500,450,375 


0.9,0.7,0.8 


NGC 4278 


2002 Feb. 


BVI 


600, 450, 360 


1.4,0.9,0.9 


NGC 4365 


2008 Apr. 


9ri 


650, 350, 300 


0.8,0.8,0.8 


NGC 5846 


2009 Apr. 


gri 


2760, 1606, 1350 


1.0,0.6,0.8 


NGC 7457 


2009 Oct. 


BVR 


6300, 6000, 7200 


0.7,0.8,0.7 



Table 2. Summary of the ground-based imaging observations. The galaxy 
name, observation date, filters employed, together with the respective 
exposure time and seeing are listed. All the observations were per- 
formed using Su baru/Suprime-Cam , except for NGC 7457 observed with 
WIYN/MiniMo iHargis etalj201ll) . 



tween NGC 1407 and NGC 1400 objects will be discussed in Sec- 
tion|5] 

Next, we calculate the radius at which the number of GC can- 
didates per unit area flattens out (see i j3.5| for the method), that is 
an estimate of the radius at which the contribution of the contami- 
nants (Glactic stars and high redshift galaxies) becomes dominant. 
Therefore, we count out of the GC selection all the objects outside 
this background radius. We have tested that the effect of the con- 
taminants on the GC colour distribution is minimal (< 0.1) mag 
and that it does not affect co nsiderably the GC colour bimodality 
(see also. lArnold et alj|201ll) . 

Finally, we perform a visual check to make sure that no out- 
liers, such as extended sources or image artefacts, contaminate the 
final GC catalogue. Most of the outliers turned out to be close to the 
galaxy (within 1 arcmin) where the galaxy light contamination and 
the crowded field makes the Subaru photometry unreliable. Within 
this radius, the contribution of the HST imaging becomes crucial. 

As far as the HST/ ACS imaging is concerned, it is worth not- 
ing that the diffraction limited quality of HST imaging has the ad- 
vantage of making extragalactic GCs partially resolved for all of 
our galaxies. Therefore, the GC selection in our space-based imag- 
ing is also based on a size selection in addition to magnitude and 
colour-colour criteria (if available). In NGC 3377 and NGC 4278, 
for which only g and z imaging are available, we flag as GC can- 
didates all the objects with colour 0.7 < (g — z) < 1.6 and sizes 
0.1 < rh < 20 pc. Such a choice is motivated by the clear drop 
off in the d ensity of the objec ts outside the adopted colour cut, as 
observed in lBlom et aT] fcOl ll) . For NGC 2768, we adopt the same 
size-cut as above and select objects within 0.4 < (B — V) < 1.2 
and 0.8 < (V — I) < 1.4 and an upper magnitude of Mi ~ —12. 
In Figure [4] the colour-magnitude diagrams of the GC candidates 
and of the spectroscopically confirmed GCs are shown. 



3.4 GC colour bimodality 

To prob e bimodality, we use a Kaye's Mixture Model algorithm 
(KMM, I Ashman et all 1 19941) that returns colour peaks, variance 
and number of objects in the detected subpopulations. KMM was 
run on all the (g — i) distributions of the Subaru GC candidates 
brighter than the turnover magnitude Mtom- In this paper we 
adopt Mxp M.i = 8.00 ± 0.03 , derived using the i — I transfor- 
m ation fromlFaifer et alJ d201 ll) to calibrate the / band TOM given 
bv lKundu & Whitmord l l200ll) into our photometric system. 
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Figure 4. GC system colour-magnitude diagrams. Photometrically selected GCs and spectroscopically confirmed GCs are shown as small and large black 
points respectively. Upper and bottom dashed lines represent the to Cen magnitude (M; ft: —11.6) and the GC turnover magnitude (AfxOM,i = — 8 mag) at 
the distances given in Table[T] 



For NGC 2768, NGC 7457 and NGC 4278, we study the 
(Rc — z), (B — V) and (g — z) colour distributio ns respectively. We 
perform a bisector fit dFeigelson & Babulll992h to the bright GCs 
in common between the HST and Subaru images, in order to con- 
vert HST magnitudes of the spectroscopically confirmed GCs into 
the respective Subaru photometric system. For NGC 4278, given 
the wide spatial coverage and better quality of its HST data, we 
transform Subaru magnitudes of the confirmed GCs into the HST 
photometric system. Results of the KMM analysis are presented in 
Table[5] We find that eight out of nine galaxies show significant GC 



colour bimodality. GCs are divided into blue and red according to 
the local minimum of the best fit Gaussians in Figure[5] 

3.5 GC spatial distribution 

We construct projected surface density profiles for the GC sub- 
populations of our galaxy sample. GC candidates brighter than the 
turnover magnitude were binned in circular annuli and corrected 
for areal incompleteness. The Subaru dataset was supplemented 
with HST photometry, in order to probe the GC surface density 
in the inner regions. An important caveat to bear in mind is that we 
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Figure 5. Colour histograms of our galaxy sample. GC candidates and spectroscopically confirmed GCs are shown as empty and grey filled histograms 
respectively. For NGC 1407, we only show GCs with % < 23. Solid lines show the sum of the KMM fits of the two subpopulations. No significant bimodality 
was detected in NGC 7457. The dearth of red GCs in NGC 3377 is due to the fact that we have excluded all the GC candidates within 1 arcmin from the 
centre. The magnitude of these objects turned out to be corrupted by the underlying diffuse stellar light. 



do not correct the Subaru dataset for magnitude incompleteness. 
This effect becomes important especially for objects fainter than 
i ~ 23.5 in the innermost regions (~ 1 arcmin) where the detec- 
tion of sources is inhibited by the su rface brightness contribution of 
the host galaxy dFaifer et alj|2oT ]]). Conversely, the completeness 
of the HST photometry in z (V) band is typically above 90 (70) 
per cent at t he turnover magnitude even in the inn ermost regions, 
respectively dLarsen et alfeOOll : I Jordan et al.ll2007l ). Therefore, we 
select all the GCs brighter than the turnover magnitude (see Fig- 
ure|4) and we use circular radial bins in common between HST and 



Subaru to correct and adjust the Subaru GC surface density to that 
of the HST data. Although the matching between HST and Subaru 
data points is arbitrary, this approach preserves the relative slopes 
of the spatial distribution of the two subpopulations and allows us 
to verify whether or not blue and red GC subpopulations show dif- 
ferent spatial distribut ions around the ho st galaxy as found from 
previous authors (e.g. jBassino et alj2006l) . 

Next, we fit the resul ting (HST + Subaru) GC density profile 
N(r) with a Sersic law dSersid [I963I) similar to that commonly 
carried out for galaxy surface brightness profiles, but in practice 
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Figure 6. Surface density profiles of the GCs brighter than the turnover magnitude. Blue and red GC subpopulations are shown as blue and red filled points 
(if from Subaru imaging) and blue and red open circles (if from HST imaging) respectively. Sersic fits or alternatively power-law fits to blue and red GC 
subpopulations are shown as blue and red lines respectively. For NGC 7457, the surface density and the fit to all GCs are shown as black points and black 
line respectively. Both the data points and the fits are background subtracted according to eq. [T] Also shown as dashed lines are the scaled and arbitrarily 
offset stellar surfa ce brightness profiles obt ained from Subaru i maging as described in Q3.2\ . For completeness, we also show the surface density profiles from 
iBlom et al.1 1201 ll) for NGC 4365 and from lHargis et al.1 feOllh for NGC 7457, respectively. 



we fit a variation of the Sersic profile from lGraham et all j2005h in 
order to quantify the background level: 



N(r) = iV e exp -6 



R 

Res 



l/n 



- 1 



+ bg (1) 



where b n = 1.9992n — 0.3271, n is the Sersic index, R c a is 
the effective radius of the GC system, N e is the surface density 
at that radius and bg is the background contamination level. In case 



a Sersic fit was not feasible due to small numbers we fit a power- 
law, Le ; _jy^r]_oc_r" that has been also used for similar analy- 
ses dSpitler etai"1l2008l) . In Figure|6]the background subtracted GC 
surface density profiles for our galaxies are shown. As found in 
other galaxies, the red GCs are more centrally concentrated than the 
blue GCs (e.gjGeisler et aljl996t iBassino et alj2006t iFaifer et al.l 



l201ll : IStraderetal.ll201ll) . 

We also compare the GC surface density to the galaxy surface 
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oalaxy ID 
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"■red 


p- value 


colour 


NGC 0821 


0.84 


1.10 


0.09 


0.15 


0.034 


(ff- 


i) 


NGC 1 400 


0.79 


1.05 


0.05 


0.11 


10 -io 




i\ 
<■) 


NGC 1407 


0.85 


0.77 


0.05 


0.1 


io- 10 


to- 


i) 


NGC 2768 


0.41 


1.10 


0.15 


0.22 


IO" 5 


(Ro- 


-*) 


NGC 3377 


0.73 


0.93 


0.05 


0.08 


IO" 6 


to- 


i) 


NGC 4278 


0.93 


1.31 


0.10 


0.16 


io- 10 


to- 


*) 


NGC 4365 


0.78 


1.04 


0.05 


0.13 


10 -io 


to- 


» 


NGC 5846 


0.74 


1.00 


0.15 


0.06 


IO -5 


to- 


>) 


NGC 7457 


1.12 


1.36 


0.13 


0.17 


0.311 


(B- 


R) 



Table 3. KMM results. For each galaxy, ^tbluc ar >d Mrcd represent the mean 
of the blue and red peak respectively, whereas <r b i uc and cr red are the Gaus- 
sian a for each peak. The last column gives the p-value, that is the confi- 
dence level with which the hypothesis of an unimodal colour distribution 
can be rejected (larger confidence for smaller p-values). The last column 
gives the colour distribution used in the analysis. NGC 7457 is the only 
galaxy without significant bimodality detected. 

brightness for each galaxy, obtained from the Subaru images using 
I RAF /Ell ipse as described in i )3.2| The surface brightness pro- 
files were corrected for the local galaxy ellipticity and then shifted 
by an arbitrary constant for comparison purposes (different for each 
galaxy). Figure[6]shows that, qualitatively, the slope of the red GC 
subpopulation agrees with that of the surface brightness of the host 
galaxy. 



4 SPECTROSCOPIC OBSERVATIONS AND DATA 
ANALYSIS 

Spectroscopic observations were performed wi th the DEep Imag- 
ing Multi-Object Spectrograph (DEIMOS, iFaber et alj 1200 
mounted on the 10 m Keck-II telescope. Galaxies were targeted 
in the period between 2006 and 201 1, employing a different num- 
ber of masks for each galaxy depending on the richness of the GC 
system. Objects selected for spectroscopic follower were chosen 
according to their likelihood of being GCs, i.e. giving priority to 
the objects selected with methods as Section [3~3l The spectroscopic 
observations are summarized in Table[4] 

The large collecting area of the 10-meter Keck primary mir- 
ror combined with the ~ 16 x 5 arcmin 2 of DEIMOS is the ideal 
combination to investigate the outskirts of galaxies where GCs are 
expected to be one of the best tracers of the total galaxy poten- 
tial. For all the galaxies, DEIMOS was set up with the 1200 1/m 
grating centred on 7800 A together with 1 arcsec wide slits, al- 
lowing coverage of the region between ~ 6500 - 8700 A with a 
resolution of A A ~ 1.5 A. Raw spectra were reduced using the 
DEIMOS/spec2d reduction pipeline provided onlin e, that produces 
calibrated and sky subtra cted spectra for each slit dNewman et al.l 
l2012l : ICooperetal.ll2012h . 

We estimate the radial velocity of the GCs by measuring the 
Doppler shift of Calcium Triplet (CaT) absorption lines that charac- 
terises the infrared part of their spectra at 8498 A, 8542 A, 8662 A, 
respectively. We measure radial velocities with IRAF/Fxcor 
that performs a cross-correlation between the Fourier transformed 
science spectrum and 13 template Galactic star Fourier trans- 
formed spectra. The template spectra were observed with the same 
DEIMOS setup used for scientific spectra and they cover a wide 
range of spectral type, luminosity and metallicity (from F to M 
type). Fxcor was configured to have all the science and template 
spectra in the same wavelength range from 8300 to 8900 A with 



Galaxy ID 


Masks 


Exp. time 


NGCs 


Nstars 


N ga l 






[hr] 








NGC 082 1 


g 


9.2 


61 


10 


4 


NGC 1 400 


4 


9 


72 


(-; 


97 


NGC 1407 


9 


20 


369 


5 


11 


NGC 2768 


5 


12.7 


109 


57 


9 


NGC 3377 


1 


8.3 


126 


16 


74 


NGC 4278 


4 


8.8 


256 


11 


33 


NGC 4365 


6 


9.0 


269 


6 


19 


NGC 5846 


6 


9.1 


195 


32 


1 


NGC 7457 


2 


1.2 


21 


11 


1 


NGC 3115 


5 


14 


190 


29 





NGC 4486 


5 


5.0 


737 


116 


27 


NGC 4494 


5 


4.6 


117 


34 


108 



Table 4. Summary of the spectroscopic observations for our twelve galax- 
ies. The table lists the galaxy name and the total number of the DEIMOS 
masks used. To take into account the different seeing conditions over differ- 
ent nights, we show the effective exposure time weig hted by the mean see - 
ing conditions during the observation time as done in ICoccatoetaIlfe009t) . 
Also shown are the total number of spectroscopically confirmed GCs (in- 
cluding marginal GCs), Galactic stars and background galaxies respectively. 
The datasets for NGC 3115 and NGC 4486 also incl ude GCs from exter- 
nal datasets as described in lArnold et all 1201 ll) and IStrader et alj fcOllh 
respectively. 

the same DEIMOS spectral resolution. The radial velocity for each 
object was estimated as the mean of the radial velocity resulting 
from the correlation with each template star. The respective errors 
were evaluated by ad ding in quadrature the default error given by 
fxcor as described in iTonrv & Davisl <fl979h to the standard devia- 
tion among the stellar templates, which is an estimate of the sys- 
tematics. 

4.1 Kinematic selection criteria 

Once the final radial velocities of the GC candidates are obtained, 
we redshift-correct the spectra and perform a visual check to verify 
that the CaT lines are real and that they lie at the expected rest 
wavelength. Our selection criteria require the presence of at least 
two "visible" CaT lines (typically the brightest 8542 A and 8662 A 
lines) and of the Ha absorption line at 6563 A, if probed in the 
redshift corrected spectra. The visual analysis is performed by at 
least two members of the team and it produces a final spectroscopic 
consensus catalogue. Spectra that show a radial velocity consistent 
with a GC, but for which it was not possible to reach a consensus, 
were flagged as "marginal" and are not included in the kinematic 
analysis. These objects have usually low signal-to-noise spectra or 
suffer from bad sky-line subtraction making the line identification 
subjective. 

Besides background galaxies that can be usually spotted by 
the emission lines in their spectra, the main outliers of our spectro- 
scopic selection are Galactic stars because they show a stellar spec- 
trum with radial velocity Vrad ~ ± 500 km s . In most cases 
GCs and Galactic stars are well-separated in velocity, but for galax- 
ies with V S ys < 1000 km s _1 Galactic stars might introduce a low 
velocity tail in the observed GC candidate velocity di stribution. We 
decide d to use a friendless algorithm introduced in iMerrett et alj 
d2003l) that flags objects deviating by more than n x a from the 
velocity distribution of their N nearest neighbours. We use n = 3 
and 10 < iV < 20 (depending on the galaxy) to exclude both 
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Figure 7. GC radial velocity distributions with galactocentric radius. Confirmed GCs, marginal GCs and Galactic stars are shown as filled points, open points 
and green boxes respectively. GCs are colour coded according to their subpopulation (blue or red) membership. Galactocentric radii have been translated into 
effective radii on the top axis. Dotted curves and dot-dashed lines are the clipping envelope as defined by the friend-less algorithm (see text) and the galaxy 
systemic velocities from Table [T] respectively. Also shown as thick dotted lines are the systemic velocities of the galaxies that might contaminate the GC 
system of the target galaxy. GCs have a mean velocity similar to the galaxy systemic velocity. 



stars and possible outliers that lie outside the 3a envelope from the 
kinematic analysis. 

In Figure [7] we show the distribution of all the spectroscopi- 
cally confirmed GCs and Galactic stars in a radius-velocity phase 
space. In this plot and hereafter, galactocentric distances are ex- 
pressed as a equivalent radius, that is defined as: 



R=JqX* + — , (2) 



where q is the axis ratio defined as the ratio of minor over the ma- 
jor axis of the galaxy (see Table [T}, X and Y are Cartesian coor- 
dinates of an object in the galaxy rest frame, with the origin at the 
galaxy centre, and with X and Y aligned along the photometric ma- 
jor axis and minor axis respectively. The final spectro-photometric 
catalogues of our galaxy sample are available on-line and they in- 
clude spectroscopically confirmed GCs, "marginal" GCs, Galactic 
stars and galaxies. 
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Figure 8. Comparison of our GC radial velocities with previous studie s. For each galaxy the root-mean-squa re of the velocity difference b etween the two 
sampl es is also shown. Literature datasets come from lBrodie et al l (12005) (black points) and lLarsen et al.l i2003l) (open circles) for NGC 4365: IChomiuk et alj 
120081) for NGC 7457; |Puzia et alj fe004l) for NGC 5846. The one-to-one line is shown as a dashed line. 



4.2 Repeated GC measurements 

We searched through the literature for GCs observed in our sur- 
veyed galaxies, finding existing datasets for three galaxies. In sum- 
m ary, we have re-ob served: (i) 6 GCs of the 26 GCs confirmed 
bv lPuzia et al1 ( l2004l) i n NGC 5846 using VL T/FORS2; (ii) 7 GCs 
of the 13 GCs found bv lChomiuk eTail d2008l) in NGC 7457 using 
Keck/LRI S ; (iii) 24 GC s of the 33 GCs resulting from th e combina- 
tion of the lBrodie et alj d2005l) and lLarsen et al. I d2003l) catalogues 
of NGC 4365, both observed with Keck/LRIS. 

In Figure [8] we compare our DEIMOS dataset with literature 
studies. There is a rough agreement between external datasets and 
our DEIMOS radial velocities, although the root-mean-square of 
the velocity difference between the two datasets is larger than what 
was found with DEIMOS repeated measurements. 

We supplement our GC catalogues with external GC radial 
velocities corrected for the mean offset between the literature and 
our DEIMOS measurements. This offset is —15km s _1 for both 
NGC 7457 and NGC 5846. I n NGC 4365 the o ffset between 
DEIMOS radial velocities and iBrodie etal I j2005l) is effecti vely 
zero, and we only correct the dataset of Larsen etal] d2003h by 
-41km s" 1 . 

We test the reliability of our spectroscopic measurements by 
comparing the radial velocity of two or three GCs observed over 
different nights. Overall, we have collected 118 repeated GC radial 
velocities and these are shown in Figure [9] One-third of the total 
repeated sample is from a NGC 2768 mask that was observed twice 
with a one night offset. We find that repeated GC measurements are 
in good agreement in all the galaxies. The root-mean-square (rms) 
of the velocity difference between two observations is SV ms = 
15 ± lions" 1 . 



4.3 Kinematic analysis 

We study the kinematic properties of our GC systems using a max- 
imum likelihood approa ch. We summarize h ere the main points of 
this analysis and refer to lFoster et alj j201ll) for the details. We di- 
vide the data in radial bins and then we fit the amplitude of the ro- 
tation (V IC ,t), the velocity dispersion (a) and the kinematic position 
angle (P^lkin) simultaneously. The bin size varies from galaxy to 
galaxy (wider for larger datasets) and it was set to have roughly the 
same number of objects per bin. For the j radial bin we minimise 
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Figure 9. Repeated GC radial velocity measurements. In the top panel we 
show the radial velocity of a GC V^. a< j 1 against the radial velocity of the 
same GC observed in a different night V r ad 2 ■ Data points are colour coded 
according their host galaxy membership listed on the right. The dotted line 
is a one-to-one line. In the bottom panel we show the difference V ra d.i — 
Vrad,2 as a function of the i band magnitude of the object. The dotted line 
shows the constant V^. a< j 2 — ^rad.i = to guide the eye. Radial velocities 
from different nights are in good agreement with each other, with an overall 
root-mean-square (rms) of 15 ± 1 km s _1 . No significant trend is seen as a 
function of GC magnitude. 
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(4) 



In eq. |4] PAi, Vohs.i and A Vobs.i are the position angle, reces- 
sion velocity and uncertainty on the recession velocity for the 
i th GC, respectively. The axis ratio of GC systems is challeng- 
ing to derive directly because of low number statistics, but it 
has been constrained in galaxies with large photometric dataset s 
dKissler-Patig et al.l 1 19971 I S trader et al.l 1201 ll. iBlom et alj 1201 ll) . 
Therefore, the kinematic axis ratio gkm of the GC system was as- 
sumed to be equal to the respective photometric axis ratio of the 
galaxy light (see Table[T}. The galaxy systemic velocity V^ ya was 
fixed to the values given in Table [Tj because the GC mean veloc- 
ity is in good agreement with the galaxy systemic velocity itself. 
Uncertainties on the fits to the GC kinematics are obtained usin g 
a bootstrapping method similar to that used by ICote et af I feooih . 
We obtain 1000 "mock" GC kinematic samples for each galaxy, by 
sampling with replacement from our measured distribution. Think 
kind of kinematic modelling tends to overestimate the rotation am- 
plitude when the kinematic position angle is a free parameter. The 
correction for this bias is described in Appendix [A] 

By way of example, the best fit to eq. [4] for the GCs in 
NGC 3377 and its kinematic dataset are shown in Figure [Tol and 
Table [5] respectively. 

The kinematic method described above was used to investigate 
the kinematic properties of our surveyed galaxies in two different 
ways: 

• First, the kinematic properties of each GC system were calcu- 
lated as a function of the galactocentric radius. This approach al- 
lows us to compare the kinematics of a given GC system to different 
kinematic probes (e.g., long-slit spectroscopy and/or PNe kinemat- 
ics). We generally use between 20 and 30 GCs per bin depending 
on the galaxy. These results are shown in Figure [TTl 

• Secondly, the confirmed GCs were sorted by their colour to study 
the GC kinematics as a function of colour. Here we use moving 
colour bins of equal width (with usually 20 GCs per bin) to investi- 
gate the effect of each GC on the final fit. This analysis is indepen- 
dent of the photometric dividing colour and it has the advantage of 
testing whether or not there is a transition in the kinematics between 
the two subpopulations. These results are shown in Figure [T2l 

To better appreciate the global kinematic properties of our GC 
systems, we also construct a 2D smoothed velocity field for our 
elongated early-type galaxies (i.e., galaxies with (b/a)ic < 0.6: 
NGC 821, NGC 2768, NGC 3377, NGC 7457). These gal axies 
are all part of the ATLAS 3D sample jKrainovic et alj|201ll) and 
are indeed good candidates to compare their moderate-fast stellar 
rotation in the innermost regions with the kinematics of their GC 
system. At every position (a;, y) on the sky we compute the local 
radial velocity ii(x, y) via interpolation of the weighted average ra- 
dial vel ocities of the N neare st neighbours, similar to the technique 
used in ICoccatoetaf] ( l2009l) : 



v(x,y) 



(5) 



with the weights being the reciprocal of the square distance d be- 
tween two neighbour GCs. The 2D velocity field was smoothed 
using a Gaussian filter kernel of variable width for each galaxy 
on a regularly spaced grid. The kernel size was arbitrarily set be- 
tween 2 and 5 kpc depending the sampling of the GC system and 
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Figure 10. Rotation of NGC 3377 GC system with position angle in dif- 
ferent radial bins. From the top to the bottom, panels show the rotation for 
all, blue and red GC subpopulations respectively. Whereas from the left to 
the right, panels show the rotation in the radial bin reported on the top left 
for each panel. The photometric major axis and the systemic velocity of 
the galaxy are both from Table [T] and represented with dashed vertical and 
dashed horizontal lines respectively. The best fit to eq.|4]and ±2cr envelope 
are shown as filled and dotted lines respectively. 



on the size of the galaxy. For the sake of compar ison, we analysed 
the SAURON data-cubes dEmsellem et~aT1l2007h of the same four 
early-type ellipticals in order to reproduce their 2D stellar velocity 
fields using our technique. These results are shown in Figure [T3l 
We also compute the root-mean-square velocity: 



1 rr. 



l 



(AV r 



rad.i 



(6) 



This quantity coincides with the velocity dispersion in eq. [4] if 
the contribution of the rotation is negligible. The u ncertainty on 
the Vr ms is estimated with the formulae provided bv lDanese et al.l 
(1980). 



5 NOTES ON INDIVIDUAL GALAXIES 

In this section we briefly discuss the kinematic results for our new 
nine GC systems that are presented in Figure [TT] Figure [121 and 
Figure [T3l 



5.1 NGC 821 



NGC 821 is an isolated E6 galaxy dde Vaucouleursl fl99ll) with 
photo metric and kinemati c signatures of an edge-on stellar 
disk dEmsellem etafl 120071 ; IProctor et at] 120091) . The dark mat- 
ter content of this galax y has been d ebated in the lit e rature . 
iRomanowskv et al] d2003t) . and recently iTeodorescu et alj d2010l) 
found that t he velocity dispersion of t he PNe decreases with radius . 
In contrast, IWeiimans etafl d2009r) , iForestell & Gebhardtl J2010h 
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ID RA Dec V ra( j <5V ra( j g 5g r 5r i Si 

[Degree] [Degree] [tar 1 ] [kins -1 ] [mag] [mag] [mag] [mag] [mag] [mag] 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) 



NGC3377. 


.GO 


161 


.94713 


14.00741 


810 


8 


21.387 





.019 


20.855 


0.013 


20.606 





,014 


NGC3377. 


.GC2 


161 


.90719 


13.96616 


724 


10 


22.295 





.032 


21.588 


0.021 


21.223 





021 


NGC3377. 


.GC3 


161 


.89561 


13.97022 


729 


13 


22.596 





,032 


21.962 


0.021 


21.671 





021 


NGC3377. 


.GC4 


161 


.92600 


14.01391 


453 


17 


22.610 





.034 


22.048 


0.023 


21.864 





,025 


NGC3377. 


.GC5 


161 


.89744 


13.97086 


681 


14 


22.789 





,037 


22.037 


0.023 


21.691 





,022 


NGC3377. 


.GC6 


161 


.93466 


13.96388 


779 


15 


22.546 





,034 


21.948 


0.023 


21.653 





,025 


NGC3377. 


.GC7 


161 


.91688 


13.96030 


617 


18 


22.919 





,046 


22.286 


0.031 


22.073 





,036 


NGC3377. 


.GC8 


161 


.88981 


13.95709 


754 


5 


20.002 





,009 


19.411 


0.006 


19.116 





,005 


NGC3377. 


.GC9 


161 


.94103 


14.03665 


804 


6 


20.150 





,009 


19.632 


0.006 


19.378 





,006 


NGC3377. 


.GC10 


161 


.87097 


13.91145 


820 


7 


20.895 





,013 


20.305 


0.008 


20.094 





,008 



NGC3377.starl 162.03633 14.02335 185 11 22.174 0.026 21.750 0.017 21.581 0.017 

NGC3377.star2 161.86723 13.90198 289 18 22.547 0.003 20.090 0.020 21.887 0.002 



NGC3377.gall 161.83237 13.92902 22.250 0.025 21.595 0.015 21.370 0.015 

NGC3377.gal2 161.95037 14.08945 24.188 0.069 24.214 0.062 24.290 0.087 



Table 5. Spectro-photometric catalogue of objects around NGC 3377. The horizontal line divides the sample into spectroscopically confirmed GCs, spectro- 
scopically confirmed Galactic stars and background galaxies. The first column gives the object ID, composed of the galaxy name and the object identification. 
Columns 2 and 3 give position in right ascension and declination (J2000), respectively. Columns 4 and 5 are the observed heliocentric radial velocities and 
uncertainties respectively. The remaining columns are the Subaru photometry in gri and respective uncertainties. The full table for all the surveyed galaxies is 
available in the online version. 



and lProctor et aljf2009h found a flat velocity dispersion for the stel- 
lar component within 100 arcsec. A ph otometric study of t he GC 
system of NGC 821 was carried out bv lSpitler et al. I d2008h using 
wide-field WIYN and HST observations. They were able to detect 
significant colour bimodality only after combining ground-based 
and space-based observations. 

Our Subaru observations suffer from moderate g band seeing 
(~ 1.2 arcsec). However, we detect significant GC colour bimodal- 
ity, although the blue and the red peaks are not clearly visible. We 
confirm that the surface density profile of the GC system extends 
up to 4 arcmin, as suggested from the WIYN imaging. 

The combination of the poor g band seeing, and the presence 
of the 10th magnitude star 2 arcmin from the galaxy centre, resulted 
in a low return rate of spectroscopically confirmed GCs. In total, we 
obtained radial velocities for 61 GCs over 6 DEIMOS masks. We 
adopt a colour split at (g — i) — 1.0. This value was used to analyse 
the kinematics of blue and red GC subpopulations separately. 

We detect significant rotation only for the blue GCs and 
for a small group of red GCs. The blue GCs are found to ro- 
tate at ~ 85kms _1 along PA^m = 85+27 deg, consistent with 
the photometr i c min or axis and in agreement with that found by 
ICoccato et al.l d2009h using PNe. Such a peculiarity is clearly vis- 
ible in Figure [T2] We will discuss this feature and its implications 
in more detail later. Interestingly, we note that the direction of the 
GC and PNe rotation coin cides with an elongate d jet/outflow struc- 
ture detected in the X-ray jPellegrini et al.l2007l) . With the adopted 
colour split, the kinematic position angle of the red GCs is gen- 
erally unconstrained, except the outer red GCs that counter rotate 
with respect to the host galaxy stars. The velocity dispersion of both 
the red and the blue GCs declines with radius with a slope similar 
to that of the PNe and host galaxy stars. 



5.2 NGC 1400 

NGC 1400 has been classi fied both as a face-o n SO l ljarrett et al .1 
120031) and as an E0 (e.g.. Ida Costa et al.|[l998l) . It is the second 
brightest galaxy in the Eridanus group after NGC 1407. We have 
assumed that NGC 1407 and NGC 1400 lie at the same distance of 
26.8 Mpc. NGC 1400 is characterised by a uniformly old stellar age 
up to ~ 1.3-Rofi dSpolaor et al ,l2008h and by an unusu ally low sys- 
temic velocity (V sys = 558 km s _1 ) for its distance. iForbes et al.l 
J2006I) studied the GC system of this galaxy using Keck/LRIS in 
imaging mode, detecting significant bimodality. 

The photometric GC selection was performed within 8 ar- 
cmin from the galaxy centre to minimise the contamination from 
NGC 1407 GC system, resulting in a clear bimodal distribution, 
with the colour separation occurring at (g — i) = 0.98. The GC 
surface density was corrected for the local NGC 1407 contribution. 
This galaxy shows the steepest red surface density profile among 
our galaxies, similar to the slope of the galaxy surface brightness. 
In contrast, the profile of the blue GC subpopulation is more radi- 
ally extended and also requires a power-law fit, because a Sersic 
function does not return a satisfactory solution. 

We present here radial velocities for 34 blue and 37 red GCs 
respectively, for a total of 71 spectroscopically confirmed GCs. De- 
spite their small angular separation (~ 10 arcmin), NGC 1407 and 
NGC 1400 also have a large peculiar velocity difference (~ 1200 
kms -1 ), assuring a reliable separation of their spectroscopically 
confirmed GCs. The red GCs mimic the rotation of the stars in the 
inner regions, whereas the rotation of the blue GCs is consistent 
with zero with a marginal signature of counter-rotation in the in- 
ner radial bin. The velocity dispersion of the red GCs is in agree- 
ment with long-slit data in the region of overlap, with a slightly 
increasing trend towards the outer regions. Conversely, the velocity 
dispersion of the blue GCs decreased with the radius. 
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Galaxy ID 


Vsys.GC 


Vrms.A 


Vrms.B 




(Vrot/o")A 


(Vrot/o")B 


(V ro t/o-) R 
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[kms" 1 ] 


[kms" 1 ] 


[kms" 1 ] 


[kms" 1 ] 
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(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 



NGC 821 
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158+^ 
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11+° 18 
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0Q+ - 27 
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NGC 1400 
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NGC 1407 
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03+ - 10 


32+ - 10 
u - JZ -o.io 


285+J0 


29 . + 12() 
zz -101 


29911 8 


NGC 2768 


1338 


165±J? 


173lj7 


160 


f 17 

-13 


n 39+0.13 


34+ 21 
u -° -0.23 


43+0-15 


97l 22 


1 21 + 32 
lzl -34 


yq+28 


NGC 3377 


685 


iool 7 


lOOlg 2 


85+ 


1(1 
7 


1Q+ 013 
u - ly -0.14 


09+ - 15 
u - uy -0.14 


n co+0.20 
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275+ 77 
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NGC 4278 
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NGC 4365 
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H+ 09 


17+ - 18 
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26+ - 11 
U - ZD -0.12 




307lf 


121+26 


NGC 5846 


1706 


238+ \l 
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f 17 

-14 


02+ 15 
u ' uz -0.12 


19+ 014 
u - ly -o.io 


12+ - 14 
u - lz -0.15 


1571?" 


259l«? 


107+ 43 
1U '-56 


NGC 7457 


847 


68±l 2 








l-6»_0.40 






32411° 







NGC 3115 710 164l| 16611? 161 ^io °- 66 ±o.09 O-^-o.ll 0.61±„;1 4 3511a A7 -ll 20 -i7 

NGC 4486 1336 32811 4 296l 28 337±1| 0.071°;°^ 0.09l°;° 7 0.29l°-j 7 105l 4 ° 92l 49 123l 39 

NGC 4494 1338 99±J| 98±® 99±lf 0.601q;1 2 0.65±°1| 0.43l°; 2 3 1741^ 167±1 6 195±H 



Table 6. GC kinematic results for our galaxy sample. Column (2) shows the systemic velocity of the host galaxy obtained from the weighted average of 
the GC radial velocities of the azimuthal-complete innermost region. The remaining columns are the overall root-mean-square velocity V r ms , the rotational 
dominance parameter (V ro t/o") and kinematic position angle PA for all, blue and red GC subpopulations respectively. 



5.3 NGC 1407 

NGC 1407 is a massive E0 gal axy at the centre of the dwarf galaxy 
dominated Eridanus A group ( iBrough et al]|2006h . It shows mod- 
erate rotation along the photom etric major axis and has a weak 
ce ntral AGN |Z hang etal.ll2007h . The stellar population analysis 
of ISpolaor et al. I J2008h found the galaxy to possess a uniformly 
old age with in ~ 0.6R c g. A dynamic al analysis of NGC 1407 
was given bv lRomanowskv et al.l J2009t) . They used Suprime-Cam 
imaging and DEIMOS spectra for 172 GCs and found a massive 
dark halo. 

In this work, we us e the photometric results presented in 
iRomanowskv et al] j2009h in which the Suprime-Cam imaging was 
reduced and analysed with the same methodology described in 
this paper. The colour distribution shows clear bimodality with the 
colour separation occurring at (q — i) = 0.98. We supplement 
the spectroscopic sample of IRomanowskv et alj d2009t) with 6 ad- 
ditional DEIMOS masks, that make NGC 1407, with a total of 
369 spectroscopically confirmed GCs, the most populous spectro- 
photometric dataset in our galaxy sample. 

We detect rotation along the photometric major axis for both 
GC subpopulations. The blue GCs rotate in the innermost radial 
bin along the photometric minor axis and between 3 and 5 R e g 
along the photometric major axis. Similarly, the rotation of the red 
GCs occurs along the photometric major axis directions, but in this 
case the rotation signal is generally larger and better constrained 
than that of the blue GCs. Both GC subpopulations appear to ro- 
tate in the outermost regions of the galaxy along the major axis of 
the galaxy. The velocity dispersion of the red GC su bpopulation in 
the inn er r egions is cons i stent with stellar data from lSpolaor et all 
j2008l) and lProctor elaH j2009h showing a decreasing profile up to 
10 i? c ff- Conversely, the velocity dispersion of the blue GCs in- 
creases with radius. 



5.4 NGC 2768 



NGC 27 68 is classified as E6 in Ide Vaucouleur sj dl99ll) and as 
SO1/2 in ISandage & B edke ( 19 94]). Spectr oscopic studies of the 
very central regions 1 McDermid et alJ uOOS) have also suggested a 
young (2.5 Gyr) stellar population associated with the disk, as sup- 



ported by the recent supernova SN 2000ds dFilippenko & Chornocld 
l200d) . iKundu & Whitmorel J200lh studied the HST/WFPC2 pho- 
tometry for 1 13 GC candidates in this galaxy, finding a statistically 
significant probability of it having a bimodal colour distribution. 

In this work, NGC 2768 was imaged with Rciz filters, 
in good seeing conditions but with the lowest exposure time 
among our sample of galaxies. Using KMM we found the colour 
distribution to be bimodal, with the blue and red peaks at 
{Rc - z) = 0.41 ±0.01 and 0.70 ±0.02 respectively. Bimodal- 
ity was also found in the HST/ACS West pointing used to design 
the DEIMOS masks in the central region. The colour peaks occur 
at(F435W-F814W)= 1.77±0.03 and2.17±0.01 respectively. 

We find 109 spectroscopically confirmed GCs over 5 
DEIMOS masks observed in sub-arcsec seeing conditions. Using 
a colour split at (Rc — z) — 0.57 and excluding the marginal GCs, 
we investigate the kinematics of the resulting 60 blue and 42 red 
GCs respectively. 

We find significant rotation only for the red GCs, which rotate 
roughly along the photometric major axis in agreement with the 
host galaxy stars. The 2-D velocity field of the red GCs is akin to 
that of ATLAS 3D in the inner regions (Figure [13). The blue GCs 
have a marginally higher velocity dispersion profile than the red 
GCs. 



5.5 NGC 3377 

NGC 3377 is a classical E6 elliptical in the Leo I group and the 
closest galaxy in our sample. T he resolved stellar properties of 
NGC 3377 have been studied by dHarris et al1l2007l) that found no 
eviden ce for any young (< 3 Gyr) stellar population. ICoccato et al .1 
j2009h studied 159 PNe in this galaxy, detecting significant rota- 
tion within 2 arcmin and a kinematic major axis twist of ~ 80 
deg fr o m the major towards the min or axis. IKundu & Whitmorel 
d200lh . lChies-Santos et al1 l l201ll) and lCho et afl < l2012h studied the 
GC system of this galaxy using HST imaging, finding the colour 
distribution to be likely bimodal. 

NGC 3377 was imaged with gri filters and gz filters for 
ground based and space based observations respectively. We found 
a high probability for the colour distribution to be bimodal, both in 
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Figure 11. GC system kinematics compared to stellar data. For each galaxy (labeled on the top right of each panel) the plots represent the tern of parameters 
(PAki n , Vrot , o") that minimise eq.|3]as a function of the radius. Data points are joined for clarity, with the error bars representing the 68 per cent confidence 
intervals. The dashed horizontal lin es in top and middle panels represent the photometric PA ± 180 deg from Table [T] and V TO t = 0, respectively. All 
the PNe data ( green stars) ar e from ICoccato et all J2009h . Integrated stellar light data are shown as small black points (if from long- slit) and large black 
points (if fr omlProctor et alj|2009i). Re ferences for long -slit data are: [Forestell & Gebhardtl feplOh for NGC 821;|Proctor et alj fe009l) for N GC 1400 and 
NGC 1407; fFried & Ilringworthl Jl994l) for NG C 2768: ICoccato et alj J2009t) for NGC 33 77; Ivan der Marel & FramJ Jl993l) for NGC 4278: bender et alj 
1 19941) fo r NGC 4365;lKronawitter et all feOOOl) for NGC 5846; ISimien & Prugniell feOOd) for NGC 7457. Black trangles in NGC 821 are SAURON-IFU 
data from lWeiimans etal]d20O9|) . For NGC 7457, because of the low number statistics, we only show the total GC kinematic profile in black. 



our Suprime-Cam and in the supplementary A CS archive imaging, 
as already found in lChies-Santos et"al] d201 lb . 

In this work, we present radial velocities for 126 GCs, ob- 
served in 4 DEIMOS masks. It is worth noting that we extend the 
kinematics of this galaxy to 8 R a s ~ 25 kpc, two times further 
than the PNe studies. We rule out two GCs that are likely to be- 
long to the spiral galaxy NGC 3377A (V sys ~ 573 km s _1 ) that 
lies ~ 8 arcmin North- West from NGC 3377. With a colour split 
at (g — i) = 0.88 we study the kinematics of 57 blue and 60 red 
GCs, respectively. 

We detect significant rotation for the red GC subpopulation 



along the photometric major a xis in agreement with PNe and long- 
slit data dCoccato et ai]|2009t) . As found for the PNe, we discover 
that the kinematic position angle of the red GC subpopulation 
twists with radius from 250 deg to 150 deg (see also Figure [Toll. 
The 2-D velocity field of the red GCs is consistent with ATLAS 3D 
in the inner regions (Figure [T3ll. The position angle of the blue GC 
subpopulation is unconstrained, implying an overall null rotation 
with radius. The velocity dispersion of the red GC subpopulation 
is very flat and it matches the stellar and the PNe. The blue GCs 
have an overall higher velocity dispersion profile than the red GC 
subpopulation. 
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Figure 12. GC kinematics as a function of colour. Kinematic position angle, rotation velocity and velocity dispersion are shown in the top, central and bottom 
panel respectively. Dotted lines represent 68 per cent confidence intervals. Also shown are the photometric major axis (horizontal dotted lines) from Table[T] 
and the colour split between blue and red GC subpopulations (vertical dashed line) as derived from the KMM analysis. 



5.6 NGC 4278 

NGC 4278 is an elliptical (El-2) member of the Coma I cloud. 
It has been extens ively studied at different wavelengths: in radio, 
iNagar et"al] d2005l) detected tw o sub-parsec je ts; in X-rays, a long 
(months) tim e scale variability dHo et alj200ll) and a dominant nu- 
clear source dYounes et al]|2010h were detected. It has long been 
known f or its massive (~ 10 8 Mq) HI disc extending beyond 
10 R eS dGallagheretaLlll977t) and used to infer the dark matter 
content at large radii l lBertola et al.ll 1993b . The bimo dality of the 
GC sy stem of NGC 4278 has been und er debate. Both lForbes et al.l 
dl996h and lKundu & Whitmorel l l200ll) investigated the GC system 
of NGC 4278 with ffiT/WFPC2, finding a marginal or "likely" 
probability for the (V — I) distribution to be bimodal. Recently, 
IChies-Santos et al.l d201lh combined HSTI 'ACS and WHT/LRIS 
imaging, finding evidence for bimodality in space-based but not 
ground-based datasets. 

In this work we present new wide field BVI Suprime-Cam 
imaging (with the B band in poor seeing conditions), combined 
with four g and z HSTI 'ACS pointings downloaded from the Hub- 
ble Legacy Archive (see Section [T2l . We find statistically signif- 
icant bimodality both in our ground-based and space-based obser- 



vations, with the colour separation occurring at (B — I) 
and (g — z) = 1.1, respectively. 



1.9 



We spectroscopically confirm 256 GCs over four DEIMOS 
masks observed in good seeing conditions. Given the wide field 
coverage of the four HST pointings, we decided to use HST magni- 
tudes to divide the kinematics of blue and red GC subpopulations, 
converting (V — I) colours into (g — z) if only Subaru photom- 
etry was available. As for other group members, we exclude from 
the kinematics analysis GCs that might be associated with other 
galaxies. We select and remove three GCs whose position and ra- 
dial velocity are consistent with the galaxy NGC 4283 (V^ ys = 
984 km s" 1 ) that lies 3.6 (~ 15 kpc) arcmin from NGC 4278. 



We find no clear evidence for red GC rotation, but there is a 
hint that the blue GCs rotate in the outer regions along a direction 
intermediate between the major and the minor axis. The velocity 
dispersion of the red and blue GCs seem to be in good agreement 
with long-slit data in the inner regions although there is no direct 
overlap between the two datasets. The blue GCs have an overall 
higher velocity dispersion at intermediate radii. 
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Figure 12 - continued 



5.7 NGC 4365 

NGC 4365 is a massive elliptical (E3) behind the Virgo cluster 
with signs of an ongoing merger (Mihos e t al. 2012, in prepara- 
tion). It has a kinematically distinct core dKrainovic et al. Him 
and an unusual ste llar rotation along the photometric minor axis 
(Ben der et al .1 19941) . Its GC system is mainly known for an odd GC 
colour distribution that seems to suggest the presence of three, in- 
stead of the classic two, subpopulati ons l lLarsen et alj|2003l . l2005l ; 
iBrodie et alj2005l : lBlom et alj|201 ll) . 



We make u se of the photometric results presented in 
iBlometalJ < l201ll) based on gz HST/ ACS and gri Suprime-Cam 
imaging. In this work, we will treat the GC system of NGC 4365 
as made up of two classic subpopulations and we refer to Blom et 
al. (2012, in preparation) for a detailed study of the kinematics of 
the three subpopulations. We use a colour split at [g — i) ~ 0.91 
that includes the "green" GCs as part of the red subpop ulation. We 
repres ent here the surface density profile published in lBlom et al.l 
d20 1 lh that shows that red GCs is more centrally concentrated than 
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NGC821-BlueGCs 



NGC 2768-Red GCs 




Figure 13. Two-dimensional smoothed velocity fields of our disky early- 
type galaxies. Black points and ellipses represent the locations of the spec- 
troscopically confirmed GCs for a given population and the 2, 4, 6 R E ff 
isophotes corrected for the galaxy ellipticity respectively. The white boxes 
are the maximum red-shifted and minimum blue-shifted velocities in the 
galaxy rest frame in unit of km s — . The SAURON stellar velocity maps 
are also shown in the innermost region of the galaxies. The agreement be- 
tween the two datasets (except for NGC 821) is extremely good. 



the blue GC subpopulation, with a trend similar to that of the galaxy 
stellar light (Figure[6]l. 

We spectroscopically confirm 269 GCs over 6 DEIMOS 
masks. We study the kinematics of 87 blue and 164 red GCs 
respectively. We a lso a dd 9 GCs from the combined dataset of 
lLarsen et al] j2003b and lBrodie et al.l d2005l) as described in ^4~2l 
We find that the red GCs rotate within R < 200 arcsec along the 
photometric minor axis, mimicking the kinematics of the stars. The 
rotation of the blue GCs is only significant at intermediate radii. 
The velocity dispersion profile is identical for the blue and the red 
GCs. 



5.8 NGC 5846 

NGC 5846 is the brightest member of a galaxy g roup. The dynam 



icsof NGC 5846 has been modelled within lf? e ff I Cappellari et al .1 



l2007l : lKronawitter et al.|[2000h . Recently. iDas et al.l d2008T > derived 
the mass of this galaxy up to 117? ff exploiting X-ray and PNe 
observations. An HS T 7WFPC2 an a lysis of the NGC 5846 GC 
system was given in iForbes et al.l i 19971) and then revisited by 
IChies-Santos et al.l feOOot) . In addition to the classic bimodality, 
they also discovered that the GC system is better aligned with the 
galaxy's minor axis than its major axis and that this galaxy has an 
unusually low specific frequency compared to similar dominant el- 
lipticals in groups or clusters. 

In this work we present new Subaru photometry in gri fil- 
ters, with the g band in moderate seeing conditions (~ 1 arcsec). 
Although the red peak is not clearly visible (Figure [5}, KMM re- 



turned a high probability for the colour distribution to be bimodal 
with a colour split at (g — i) — 0.95. 

We spectroscopically confirm 195 GCs over 6 DEIMOS 
masks. We also note that some radial velocities might be asso- 
ciated with other bright members of the group as NGC 5846A 
(V sys = 2200 km s" 1 ) and NGC 5845 (V sys = 1450 km s" 1 ). 
However the similar systemic velocities of these galaxies with 
NGC 5846 (V sys = 1712 km s -1 ) makes any attempt to distin- 
guish their GC populations problematic. We also sup plement our 
DEIM OS catalogue with 22 GCs from the dataset of IPuzia et al .1 
I l2004|) as described in i j4.2l In summary, we study the kinematics 
of 104 blue and 108 red GCs. 

We detect rotation only for the blue GCs, which rotate be- 
tween 150 and 300 arcsec close to photometric major axis. This 
feature might be caused by the GCs of NGC 5846A, which con- 
taminates our sample within 300 arcsec. The velocity dispersion of 
the red GCs is flat with radius and consistent with other studies, 
although the PNe seem to suggest a slightly decreasing slope. The 
blue GCs have a systematically higher velocity dispersion than the 
red GC subpopulation. 

5.9 NGC 7457 

NGC 7457 is an isolated SO with a pseudo-bulge dete cted both pho- 
tomet rically dTomita et al.ll2000T) and kinematically jPinknev et al.l 
120031) that sho ws an unusually l ow cen tral velocity dispersion for 
its luminosity. lEmsellem et al.l J2007t) revealed that this galaxy 
has a s mall counter rotating core that might be the result of a 
merger. Ichomiuk et al.l J2008b gave an overview of the GC system 
of NGC 7457 using HST observations and Keck/LRIS spectra for 
13 GCs. They find evidence for a third intermediate population of 
GCs sharing the same age (2 -2.5 Gyr) as the youn g nuclear (radius 
of 1.5 arcsec) stellar population discovered by Sirchenko et al .1 



NGC 7457 is the only galaxy in our sample not observed with 
Suprime-Cam. We have instead used the photometric results of 
lHargis et all d201lh to design our DEIMOS masks, and we refer 
to their paper for a detailed description of the data reduction and 
analysis. They provide wide field WIYN/Minimosaic photometry 
of the GC system of this galaxy, finding that the colour distribution 
is not bimodal and showing the total population to be made up of 
210 ± 30 GCs with a radial extent of 12 ± 2 kpc. 

We find 21 spectroscopically confirmed GCs in two DEIMOS 
masks observed in average 0.9 arcsec seeing conditions. We 
also obtained spectra for 7 of the 13 GCs already observed by 
IChomiuk et all 00081) . finding their LRIS and our DEIMOS radial 
velocities to be consistent within the errors, but with a mean differ- 
ence of ~ 12 km s _1 (see £|4.2| l. Therefore, we add the remaining 
6 Chomiuk GCs corrected for the velocity offset to our GC sample. 
The final GC catalogue consists of 27 GC radial velocities. 

Our confirmed GCs appear to be distributed within 2 ar- 
cmin ~ 7.7 kpc from the galaxy centre, reflecting the expected 
poor extent of the underlying GC subpopulation. The mean veloc- 
ity of the system (V^ YS = 847 ± 20 km s -1 ) is in good agreement 
with both lChomiuk et al.l (2008) and literature data. As this galaxy 
appears to be unimodal in colour, we fit the kinematics of all the 
GCs without any s plit in colour and we com pare the results with 
the long-slit data of lSimien & Prugniell l l2000h . We find that the GC 
system of NGC 7457 is rapidly rotating along the photometric ma- 
jor axis, with a flat rotation curve at ~ 80 km s _1 up to 2R e ff. In 
contrast, the velocity dispersion is low (< 50 km s _1 ) at all radii. 
The agreement with the ATLAS 2D velocity map is also very good 
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(Figure[T3V We will discuss the implication of these results in more 
detail in fll04l 



6 GC FORMATION MODELS 

Having determined the photometric and kinematic properties of our 
GC sample, in the following sections we will discuss these results 
as a whole. We look for common kinematic features that might re- 
tain key information about the formation history of the galaxies 
themselves and then we compare these results to theoretical mod- 
els. In order to do so, it is first worth giving a brief summary of the 
main GC formation scenarios proposed in the literature, focusing 
on their kinematic predictions. 

To date, the formation of GC bi modality has been investi- 
gated both in a cosmological context jKravtsov & Gnedin|[20o3 ; 
iDiemand et al.l2005l ; lMoore et al.l2006l) and at smaller scales, with 
models fine tuned to reprod uce the properties of specific galax- 
ies (see lVesperini et alj|2003l for NGC 4486 or lDeason et ailhoill 
for the Milky Way and M31). The challenge of producing simu- 
lations of GC bimodality formation is intertwined with our poor 
understanding of how GC systems formed in the first place (e.g., 
iBeaslev et alj|2002l ; lElmegreen et al.ll2012h . As a consequence, the 
three classic formation scenarios proposed in the literature and 
summarized below have few, if any, theoretical predictions directly 
comparable with the observations. In particular, there is a general 
dearth of GC kinematics predictions. 

In the major merger scenario I Ashman & Zepjl 19921) two or 
more gas-rich disk galaxies with pre-existing blue and red GC sub- 
populations merge to form an elliptical galaxy. New r ed GCs may 
form from the star formation induced by the merger iBekki et"aU 
l2002h . 

This scenario has been simulated by |Bekkietai] ( l2005l) (here- 
after B+05). In this simulation the pre-existing blue and red GCs 
are assumed to have a Milky Way-like spatial distribution, but they 
are both pressure-supported. This is a reasonable assumption for 
the blue but not for the red GCs of the Milky Way that are known to 
be rotation-supported dCotell 19991 ; iDeason et alJbo 1 ll) . B+05 pre- 
dicted that a merger with even mass-to-mass ratio produces strong 
rotation for both GC subpopualtions in the outer regions of the rem- 
nant. This result is independent of the orbital configuration of the 
merger. The velocity dispersion profile is predicted to decrease with 
radius, but it would be flatter in case of multiple mergers. Also, the 
blue GCs are expected to show a larger central velocity dispersion 
than red GC subpopulation. Finally, the ratio of the maximum ro- 
tational velocity V m to the central velocity dispersion ao of the GC 
systems ranges from 0.2 to 1.0 within 6R c a for both the blue and 
red GC subpopulations, but mo st GC systems, viewe d from various 
angles, have (Kn/oo) < 0.5. iBekki & Peng) d2006l) carried out a 
simulation similar to that of B+05 (with the rotation of the disk in 
the spiral progenitors included), fine tuned to study the dynamics 
of the planetary nebulae (PNe) in elliptical galaxies. They found 
that the effect of the residual spin disk (additional to the initial or- 
bital angular momentum) on the final PNe kinematics enhances the 
rotation at all radii, making the merger remnant rotation-supported. 

In the mu ltiphase dissipational collapse scenario 
dForbes et ail 1 19971) . GCs are the result of an early two phase 
collapse, and hence two main star formation episodes, that the 
galaxy undergoes. The blue GCs form during the first star forma- 
tion episode in a metal-poor cloud, whereas red GC subpopulation 
form in a second phase after the gas in the galaxy is self-enriched. 
In this scenario, a fraction of blue GCs might also come from 



accreted satellite galaxie s, similar to the most recent two-phase 
galaxy formation model l lOser et al.|[201oh . No numerical simula- 
tions have been performed so far to test this model. Since the red 
GCs a re coeval with the galaxy stellar component, iForbes et al.l 
i 19971) infer that they should share the same kinematic properties 
as the galaxy stars, and hence also the same spatial distribution. 
No significant rotation is expected for the blue GC subpopulation . 

In the dissipationless accretion scenario ( tCote et alj|l998h . 
the red GCs form from the monolithic collapse of a giant proto- 
galactic cloud, whereas the blue GCs are accreted from low mass 
galaxies due to mergers or tidal stripping. To date, s imulations 
based on similar scenarios have been performed (e.g., lOser et all 
l20ld) but they do not include GC particles, and they have not made 
kinematic predicti ons. In add i tion to the observation of infalling 
co-moving groups, ICote et alj fl998h also infer that we should ex- 
pect to observe blue GCs with radially biased orbits without overall 
rotat ion. 

IBekki et all J2008l) (hereafter B+08) performed a high resolu- 
tion N-body cosmological simulation combined with semi-analytic 
models. They found that almost half of the simulated 10 5 galaxies 
show clear bimodality in their GC metallicity distribution. The ma- 
jority of GCs form in low-mass galaxies at redshift greater than 3 
with the blue GCs being slightly older (< 1 Gyr) than the red GC 
subpopulation. B+08 also made predictions for the overall kine- 
matics of each GC system, whereas current spectroscopic observa- 
tions only sample a smaller fraction of it. They find that the velocity 
dispersion of blue and red GC subpopulations increases with the to- 
tal luminosity of the host galaxy and that the ratio of the velocity 
dispersion of the blue to red GC subpopulation is (ob/or) ~ 1 
for a wide range of luminosity. Finally, they also predicted that the 
GC systems of most galaxies are mainly pressure supported with 
(Vmt/c) < 0.3. This is in contrast with the wider range of Kot/ff 
found by B+05. 



7 GLOBAL RESULTS FROM OUR GC SAMPLE 

In this section we summarize the results for the GC systems pre- 
sented in this paper, including the three galaxies previously pub- 
lished by us (i.e., NGC 3115, NGC 4486, NGC 4494). Section [8] 
will incorporate other literature data, and tackle topics not covered 
here. 

7.1 Spatial Distribution 

Significant colour bimodality was detected in all the GC systems, 
except for NGC 7457. All GC system formation scenarios, as 
well as numerical simulations, predict that this bimodality should 
be reflected in different spatial distributions of the two subpop- 
ulations around the parent galaxy. A glance at Figure [6] reveals 
that this is also the case for our galaxy sample. Generally speak- 
ing, the red GCs are more centrally conc entrated than the blue 
GCs, as already found in other studie s (e.g. JSchuberth et"al]l2010l ; 
iFaifer et alj|201 ll ; IStrader et al .11201 ll) . 

The advent of wide field imaging has shown that the radial 
distribution of the red GC subpopulation matches that of the host 
galaxy light, su ggesting that they m ight have shared a similar for- 
mation history fBassino et al J 1 20061). Here we have exploited the 
34 x 22 arcmin 2 field of view of our Suprime-Cam imaging to 
compare the GC surface density to the surface brightness of the re- 
spective host galaxy (shifted by an arbitrary constant). Qualitatively 
speaking, the galaxy starlight has a similar slope to that of the red 
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Figure 14. GC rotation velocities and velocity dispersions compared to the galaxy stars/PNe. Rotation velocities (left panel) and velocity dispersions (right 
panel) are extracted within 3i? c ff . The exact extraction radius is usually set to 2R c ff but this threshold is arbitrarily adjusted by ±1-R c ff in case the GCs 
do not overlap with the stellar data. The blue and the red GC subpopulations are shown in the bottom and upper panel respectively. Dashed lines represent 
a one-to-one relationship. Blue GC data for NGC 4365 and NGC 4278 are not shown due to limited radial range. NGC 7457 is not shown. The red GC 
subpopulation shows better agreement with stars than the blue GCs. 
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Figure 15. Velocity dispersion profiles. For each galaxy, we show the ra- 
tio of the velocity dispersion of blue to red GC subpopulations normalised 
to 3 -R c ff (vertical dashed line) as a function of the effective radius. Dif- 
ferent colours correspond to different galaxies. The horizontal and vertical 
lines show (tb/cr) = (ctb/cr^ /J eff and R = 3iJ ft, respectively. 
NGC 7457 and NGC 4494 are not shown. 



GC subpopulation. The discrepancies between the surface bright- 
ness and the GC spatial distribution in innermost regions (e.g., see 
NGC 4365 and NGC 1407 in Figure|6j is due to the core-like dis- 
tribution of the GCs, that makes the GC density profile flatter. This 
feature is interpret ed as the effect of the GC disruption, stronger i n 
the central regions dBaumgardtl 19981 ; lBaumgardT & Makinol 2003h . 



7.2 Rotation 

We find a large variety of GC rotation profiles. Both the blue and 
red GCs show some degree of rotation, but there does not seem to 
be a clear correlation between the rotation patterns and the property 
of the host galaxy. A glance at Figure [TTI reveals that the red GCs 
rotate more consistently with the photometric major axis than the 
blue GC subpopulation does. For our galaxy sample the photomet- 
ric position angle coincides with the photometric major axis of the 
galaxy stars (except for NGC 4365), which means that the rotation 
velocity of the red GC subpopulation is similar to that of the galaxy 
stars. This is also true for the overall GC system of NGC 7457. 

To quantify this phenomenon, in Figure Q4] we compare the 
major axis rotation velocity of the blue and the red GC subpopu- 
lations within 3 -Rett to the rotation velocity of the stars and PNe 
(if available) at the same galactocentric distance. The exact extrac- 
tion radius was set to 2R c g, but in some cases this was relaxed by 
±li? c ff to maximise the overlap between the GCs and the stars. The 
rotation velocity of the host galaxy stars was extracted along the 
photometric major axis for all the galaxies, except for NGC 4365 
in which the bulk of the stellar rotation is occurring along the pho- 
tometric minor axis. The red GC rotation for our galaxies is similar 
to those of the host galaxy stars, at least for galaxies with a con- 
spicuous amount of stellar rotation. On the other hand, the star-GC 
system coupling is also evident for only blue GC subpopulations, 
such as in NGC 4494, NGC 821 and NGC 3115. 

The rotation of the blue GCs is more puzzling than the red 
GCs. Overall the rotation of the blue GCs is lower than the red 
GCs and often consistent with zero, yet we also detect minor axis 
rotation in NGC 821, also seen in PNe, in contrast with integrated 
starlight. 

We detect GC rotation at large galactocentric radii for some 
galaxies. This feature is only significant in NGC 1407, NGC 4486 
and only marginal in NGC 4278. The outer rotation occurs in the 
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Figure 16. Summary of GC system kinematics as a function of the host galaxy properties. Literature data and our sample are shown as red and black 
points respectively. From the top to the bottom: root-mean-square velocity for the overall GC system (V rm s a)> ra ti° of the V rms of the blue to the red 
GC subpopulation (Vrms.B/Vrms.R.), ratio of the rotation ve locity to the velocity dispersion for the overall GC system (Vrot/c)A- From the left to the 
right: X-ray luminosity (l og Lx) from lO'Sullivan et alj feOOlh at the distance given in Table[T]and Table|7] central stellar velocity dispersion (<ro iS tars) f rom 
IPrugniel & Simienl fl99^) and total K-band magnitude (M K ) from 2MASS. The black and red arrows indicate the position of NGC 7457 at (V ro t /a) A = 1-68 
and of NGC 3311 at V Tias a = 653 km s -1 , respectively. The horizontal dashe d line marks K- ms . B/Vrms,R = 1 to guide the eye. The grey shaded shaded 
region represents the theoretical predictions from the cosmological simulation of lBekki et al l bOOa) . 



same direction for both the blue and the red GC subpopulations and 
it usually coincides with the photometric major axis. 



7.3 Velocity Dispersion 

Our results show that the velocity dispersion profiles are quite flat 
for most of the galaxies, both as a function of radius (FigurefTTTl and 
as a function of the colour (Figure [T2l, The slope of the velocity 
dispersion profiles (including the contribution of the rotation) will 
be discussed below. 

As for the rotation, we study how the velocity dispersion of 
the two GC subpopulations compares to the velocity dispersion of 
the host galaxy stars. We extract the velocity dispersion at the same 
galactocentric distances as done with the rotation. Figure [T4l shows 
that the velocity dispersion of the red GC is very similar to that of 
the host galaxy stars in the region of overlap. Conversely, the blue 
GC subpopulations seem to avoid the one-to-one line. The standard 
deviation from the one-to-one line is 22 km s _1 and 50 km s _1 for 
the red and the blue GC subpopulations respectively. 

To quantify the differences between the velocity dispersion 



profiles of the two GC subpopulations, in Figure [T5] we plot the 
ratio of the velocity dispersion of the blue GCs to the velocity dis- 
persion of the red GCs (<tb /ctr) normalised to 3 R c ft as a function 
of the effective radius. In this case we show the velocity dispersion 
using rolling radial fits to appreciate the details of the profiles. Fig- 
ure 1 151 shows that the ratio ratio (o"b/o"r)/(ctb/o"r)3^ is gen- 
erally consistent with 1, but it increases towards the outer regions 
for some galaxi es, as already observed for more massive galaxies 
dLee et alfeoich . 

The overall velocity dispersion profiles are in rough agreement 
with their merger remnants assembled via multiple mergers, with 
the velocity dispersion of the blue GCs being larger or equal to that 
of the red GCs. However, the model fails in reproducing the irregu- 
larity of some dispersion profiles (e.g., NGC 1407 and NGC 4278). 



8 INCLUDING LITERATURE DATA 

We examined the literature for galaxies with a reasonable number 
of GC radial velocity measurements. We require at least 20 spec- 
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Table 7. Kinematical properties of literature GC systems. The host galaxy name (1) and Hubble Type (2) are from NED database. K band absolute magnitude 
(3) is from 2MASS apparent magnitude at the distances given in column 4 and corrected for the foreground Gal actic extinct i on fro m NED database. The 
distance in Megaparsec (4) was obtained by subtracting 0.06 mag ( Mei et al. 2 0071) from the distance modulus from lTonrv et all i200lh . except for NGC 33 1 1 
for which we adopt the distance of the Hydra I cluster of (m — M) = 33.37 mag iMieske et al.l2003) . The number of spectroscopically confirmed GC system 
(5) is after the outliers filtering as described in the text. Columns (6), (7) and (8) are the root-mean-square velocity V rms for all, blue and red GC sub populations 
respe ctively. Columns (9), (10) and (1 1) are t he (V/cr) value for a ll, blue and red GC subpop ulations respectively. The dat a referen ces are: P+04 jPuzia et alj 
12004): S+10 iSchuberth et al.l20ld): M+ll <Misgeld et alj201ll): P+ 06 IPierce et all2006l): B+06 <Bergond et alj2006l) : N+12 iNorris et alj2012h : C+03 
ICote et al.l2003l) ; H+08 iHwang et alj2008h ; L+08 iLee et al.l20olj) ; W+10 iWoodlev et al.ll201fj) . 
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Figure 17. GC system V Tnls as a function of the host galaxy properties. 
Literature data and our sample are shown as red and black points, respec- 
tively. The top and the bottom panels represent the velocity dispersion of the 
red and blue GC subpopulations as a function of the X-ray luminosity (log 
Lx), central stellar velocity dispersion (cro lS tars) and total i\"-band magni- 
tude (Mk), respectively. The solid and dashed lines represent the weighted 
linear fit to the data and its rms in km s~ 1 quoted on the bottom right of 
each panel, respectively. The dependence of red GC subpopulation on the 
galaxy properties is tighter than for the blue GCs. 



troscopically confirmed GCs for a galaxy to be included in the liter- 
ature sample. We use the same friendless algorithm as described in 
^4. 1 I to clip outliers in radius-velocity space. This selection returned 
ten GC systems, which are summarized in Table [7] Mos t of these 
galaxi es are dominant group/cluster galaxies discussed in lLee et all 
whereas the less massive galaxies are part of the survey of 



IPuzia et alj d2004h . We also include NGC 3923 from lNorris et al.1 
(2012) which shows strong signatures a recent interaction. A first 
comparison between our GC dataset and external datasets was al- 
ready given in Figure[2]in which we showed that we cover a wider 
range of mass with three times better velocity accuracy than previ- 
ous studies. 

We ran the kinematic likelihood code on the literature galax- 
ies to get the overall best fit values to eq. [4] but we do not show 
the kinematic radial profiles for the literature sample. We investi- 
gate the kinematics of blue and red GC subpopulations for seven 
out of ten galaxies using the dividing colour quoted in the re- 
spective papers. For consistency with the kinematic analysis of 
our galaxies, we do not exclude UCD candidates from the litera- 
ture data. We find that our likelihood-method tends to reduce both 
the velocity dispersion (by ~ 10 km s _1 ) and the rotation am- 
plitude (up to ~ 30 km s~ ) if the contribution of the rotation 
is low (Kot/cr < 0.4). We attribute this discrepancy to the fact 
that some previous GC kinematics studies have employed error- 
weighted least-square rotation fittings that are less appropriate for 
a system with an intrinsic dispersion. Another explanation might be 
that, to our knowledge, previous studies have never explicitly taken 
into account the bias introduced by an unconstrained position an- 
gle (see j)4.3b . which has an important effect at lower velocities. 
It is important to note that the maximum-likelihood approach does 
not alter the general results of previous works. 

A final caveat to bear in mind concerns the different kine- 
matic axes convention. Based on the fit of the GC radial velocities 
with the position angle (Figure [Toll, we define the kinematic po- 
sition angle as the angle between the direction of maximum rota- 
tion amplitude and North, that is 90 degrees away from the angular 
momentum vector. Some previous studies have instead defined the 
kinematic position angle as the direction around which the rotation 
is occurring. In other words, our kinematic maj or axis would cor - 
respond to the kinematic minor axes quoted by iLee et all d2010h . 
Our convention is in line with that used in other galaxy kinematic 
studies (e.g. ATLAS 3D ). For the galaxy properties of the literature 
sample we use the same sources quoted in Table Q] 
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Figure 18. Rotational dominance parameter of the stellar component A p 

-"■erf 

versus the galaxy ellipticity ejj . Both quantities are measured within 
1 -R ff . Data points are represented by the NGC number of our galaxy sam- 
ple (bl ack) and literature data (red). Measurements are from lEmsellem et al] 
1201 ll) if the galaxy is in ATLAS 3D .Ifnotthey are estimated from literature 
data (see text). The magenta curve is the inferred edge-on average relation 
for nearby fast-rotators. The solid green curve is (0.31±0.01) X ~ 
with the dashed lines representing its uncertainties. Galaxies are divided 
into fast-rotators and slow-rotators if they are above or below the green 
curve respectively. 



9 RESULTS 

In this section we analyze the extended sample of 22 GC systems 
(literature plus our own galaxies) , 18 of which have kinematics for 
both GC subpopulations. We carry out a number of tests to compare 
the kinematic properties of the two subpopulatinos with each other 
and with host galaxy properties. 



9.1 Correlations with host galaxy properties 

In Figure[T6]we show the GC system kinematics as a function of the 
general properties of the host galaxies. We n ow have the chan ce to 
better constrain some relations examined by iLee et al] d2010h that 
were uncertain because of their limited sample size of six galaxies 
(see their figure 12). 

The upper panels of Figure [T6] show that, as found from pre- 
vious studies, a correlation exists between Kma of the GC systems 
and their respective host galaxy X-ray luminosity Lx, central ve- 
locity dispersion of the stars <to and absolute A'-band magnitude 
Mk ■ We find that these correlations also hold for less massive 
galaxies. The correlations with galaxy properties are always tighter 
for our GC dataset as consequence of the improved reliability of 
our velocity measurements. 

The correlation of V m s with galaxy properties holds when we 
plot the Vrms of the blue and red GC subpopulations separately (see 
FigurefTTt. The main difference here is that the two subpopulations 
scatter differently with respect to the best fit lines. The rms of this 
difference for the blue GCs is systematically larger for the red GCs. 
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Figure 19. The rotational dominance parameter of GC systems. The plot 
compares the ( V TO t / c) parameter of the blue GCs (V ro t /c)b to that of the 
red GCs ( V TO t /"Or ■ The dashed line marks a one-to-one line. GC systems 
are colour coded according to the kinematics of the host galaxy: purple if 
slow-rotators and orange if fast-rotators. Filled and open circles represent 
our GC systems and the data from the literature, respectively. 



In Figure[l6] the middle panel shows that most of the GC sys- 
tems have Kms.B /Kms.R > 1. This is seen more clearly as a func- 
tion of the X-ray luminosity which is larger in galaxies with an 
extended massive halo to whi ch most of the blue GCs are thought 
to belong dForbesetal.ll2012l) . These results are in general agree- 
ment with the cosmological simulation of B+08 that predicted the 
Kms,B/Vrms,R ratio to be slightly > 1 for a wide range of host 
galaxy magnitudes. We also tested the significance of the apparent 
correlation between Kms,B/Kms,R and Mk finding a Spearman's 
rank correlation coefficient of r s = 0.35, implying that a correla- 
tion is unlikely. 

In the bottom panel of Figure [T6] we quantify the GC system 
rotation by studying the overall (Kot/c), which estimates the con- 
tribution of the GC rotation over its random motions. We find that 
galaxies are spread between < (V ro t/o")A < 0.8 except for 
NGC 7457 with (V TO t/a)A = 1-68. It also appears that galax- 
ies with lower Lx and lower oo are more rotation-supported than 
other galaxies. The bulk of GC systems with (Kot /cr) > 0.3 is con- 
sistent with cosmological simulation of B+08. However, as previ- 
ously discussed, trends of rotation with radius and colour are also 
important to examine. 



9.2 Rotation 

To make things clearer we divide the host galaxies into two kine- 
matic groups that will be used throughout this section. After con- 
sidering various possibilities, we finally decided to divide the GC 
systems according to the rotation dominance parameter of their 
host galaxy. We use the Ac> parameter within l_R e ff as de- 
fined in ICappellari et all ( 120071) with a threshold of 0.31 x y/e 
to discriminate between slow and fast-rotators, where ep „ = 
1 — yb/a) is the galaxy ellipticity. For our sample, we use 



© 2012 RAS, MNRAS 000,[TJ-?? 



Globular Cluster Kinematics in Twelve Early-type Galaxies 25 




J I I I I I I I I I I I I I I I I I I I I I I I I I I I I L 

2 4 6 8 10 12 14 

R / R eff 



Figure 20. (V^ot/c) radial profiles for our galaxy sample. Colours are as in Figure [l9l The top and bottom pa nels show the rotational dominance parameter 
for the red and blue GC systems respectively. The black lines and arrows are predictions fromlBekki et al.l (2003) an d represent the range of ( Vmax /co ) within 
2 and 6 R c (f respectively. There is no clear trend of ( V TO t / c) with radius as predicted from B+05. 



from lEmsellem et all < |201 lh if the galaxy is in ATLAS 3D , whereas 
we use the conversion formulae provided by the same authors to 
convert Vmax/cro into Ap „, where Vmax and a a are the maximum 
rotation amplitude and the central velocity dispersion of the host 
galaxy respectively (Figure [T8t. This separation has the advantage 
of being independent of the galaxy inclination. It also allows us 
to test if the slow-fast rotator separation in the inner regions per- 
sists in the outer regions. The first group of galaxies is composed 
of slow-rotator galaxies and includes mostly round massive galax- 
ies at the centre of groups or clusters. The second group includes 
fast-rotators but it covers a wide range of galaxy masses and mor- 
phologies, from E0 to SO. Fast-rotator galaxies with \p> > 0.25 
have stellar disky-like structures (photometric and/or kinematic) in 
the inner regions. 

In Figure[l9]we compare the overall rotational dominance pa- 
rameter (Vrot/c) of the blue and red GC subpopulations colour- 
coded according to their fast/slow rotator classification. We find 
that a dichotomy exists between the GC systems in slow and fast 
rotators. Most of the GC systems in the slow-rotator galaxies have 
< (Vrot/c) < 0.3. Again, we want to emphasise that this does 
not necessarily imply that these systems lack rotation because if 
their rotation occurs only in a limited radial range, it may have 
been smeared out in the overall (Kot/o"). The GC systems in fast- 
rotator galaxies are more rotation-supported and in general at least 
one GC subpopulation has a conspicuous amount of rotation. There 
is a group of galaxies with (Kot/o") ~ 0.6 which includes the fast- 
rotator galaxies with large Ajr, but it also includes NGC 4649 
with A d „ ~ 0.12. 

J >-crr 

In Figure|20]we compare the radial (Vrot/f) profiles for the 
blue and red GC systems with the numerical simulations of B+05. 
We note that B+05 quote the ratio of the maximum rotation ve- 
locity to the central velocity dispersion (Vmax/co) within 2 and 
6 Red, whereas we plot the overall (Kot/cr) in each radial bin. 



Therefore, the predictions plotted in Figure [20] are typically up- 
per limits. Another caveat is that B+05 simulate dry mergers, 
which leads to a systematic mismatch with the observed rotation 
in the galaxy centers . However, this effect is small for R > 2R c r 
dHoffman et al.l2010h . which is the radial range with better GC cov- 
erage and the most relevant for this study. 

We find that for galaxies in the fast-rotator group, at least 
one GC subpopulation has (V to t/o~) 7^ 0. Slow-rotator galaxies 
have generally slower rotating GC systems. Although the compar- 
ison with B+05 suggests an overall agreement of our results with 
their predictions, the scatter in both (V TO t/a) and in the simula- 
tions is too large to draw any strong conclusions. In more detail, 
the simulations predict (Kot/o") increasing monotonically with ra- 
dius but we see no evidence for this feature in our data. More- 
over, a direct comparison with B+05 simulations is not possible 
because they did not publish the (V I0 t/o~) profiles to 6R c g. Fig- 
ure [20] also shows that many GC systems are characterized by 
significant rotation spikes. These features might be partly caused 
by spatial incompleteness, outliers or projection effects. Alterna- 
tively, one speculative possibility is that these rotation spikes are 
the imprints of minor mergers. These events are expected (e.g., 
I Vitvitska et al.l [20021; iFfetzn ecker & Burkert 20061), and observed 
(e.g.. lRomanowskv et alj|201lh , to perturb the halo kinematics of 
massive galaxies. Detailed studies of GC syst em phase-space for 
each galaxy are needed to quantify this effect. iNorris et al.l J2012h 
also suggested that the missing outer rotation could be explained 
if the reservoir of angular momentum was located beyond the 
radii mapped by current GC system studies. This idea would be 
supported by the large rotation (T4ot/c ~ 1-0) we observe at 
~ 10-Rcff for both GC subpopulations of NGC 1407. However, at 
such galactocentric distances, results are uncertain because of spa- 
tial incompleteness and low number statistics. The same concerns 
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Figure 21. GC system position angle misalignments. The horizontal and 
the vertical axes show the misalignment between the photometric and the 
kinematic major axes for the blue and the red GC subpopulations respec- 
tively. Dashed lines mark a constant 90 degree angle. Rotation along the 
photometric major axis occurs at deg, rotation along the minor axis at 
90 deg and counter rotation at 1 80 deg. Galaxies with both a blue and red 
GC kinematic position angle unconstrained are not included for clarity. The 
rotation of the red GCs is more consistent with the photometric major axis. 



were addressed bv lStrader et al.l J201 ll) on the suspicious outer ro- 
tation (Vrot/tr ~ 0.6) observed in NGC 4486. 



9.3 Kinematic misalignment 

It is also important to study how the rotation of the GC systems 
aligns with respect to the host galaxy isophotes. An eventual posi- 
tion angle misalign ment might contain inf ormation about the triax- 
iality of the system dKrainovic et alj|201 lh . 

In Figure [21] we compare the PA misalignment of the blue 
and red GC subpopulations. For clarity, we have removed from the 
plot galaxies for which neither the PA^ n of the blue nor of the 
red GC subpopulation was constrained. A caveat here is that we 
study overall PAkm that will be biased by kinematic twists with 
radius. Broadly speaking, we note that most of the galaxies are lo- 
cated in the bottom-left quadrant of Figure[2T]towards | APA\ = 0, 
suggesting that the red GC subpopulation rotates more consistently 
with the photometric major axis than the blue GC subpopulation. 
We do not see a sharp separation between the GC systems in fast 
and slow-rotator galaxies. The GC rotation in fast-rotator galaxies 
seems to be equally consistent with the photometric major axis for 
both GC subpopulations. 

The B+05 simulations also predicted kinematic misalign- 
ments with respect to the galaxy isophotes, albeit with an amplitude 
smaller than our findings. Depending on the orbital configuration, 
these authors can also reproduce the minor-axis rotation for the blue 
and red GC subpopulations, but they do not quantify this effect. Un- 
fortunately, B+05 do not provide radial kinematic misalignments 
profiles, but they find that the rotation of both GC subpopulations 
in the outer regions of a merger remnant should be aligned along 




Figure 22. Log-slopes of the GC system V IIns profiles. Symbols are as in 
Figure [T9l For each GC system the plot shows the best fit 7 parameter to eq. 
[8]for all (top panel), blue (middle panel) and red GC subpopulation (bottom 
panel) as a function of the absolute magnitude Mx of the respective host 
galaxy. Lines are weighted linear fits to the data, whereas the dashed lines 
show an ideal flat Vrms profile (7 = 0) . The blue GC subpopulation tends 
to have a decreasing Vrms slope (7 < 0) in more massive galaxies, while 
the red GC subpopulation has steady 7 > over the whole magnitude 
range. 



the photometric major axis, which is in rough agreement with our 
findings. 



9.4 The slope of the K ms profile 

The V Tms , sometimes freely named as velocity dispersion if the 
contribution of the rotation is low, is related the line-of-sight kinetic 
energy of a galaxy per unit of mass. Its profile has been classically 
used to model dark matter halos. Here, instead of showing the de- 
tailed Kms profile for our GC system sample, we think it is more 
convenient to parametrize these profiles and to study their overall 
slopes. We use a maximum-likelihood approach to fit a power-law 
function to the non-binned Kms profile, similarly to that done in 
eq.|4]for the GC system kinematics. I n this case the \ 2 function to 



be minimized is dBergond et al.| [2006 ) 



= E 



(Vre 



+ (AKad.O 2 



+ ln(K 2 ms + (AKad, ! ) 2 ) 



and the function to be modelled is: 



V? ms (R) = VS»,o x (R/Ro)' 



(7) 



(8) 



where Ro is a scale parameter, set here to the median galactocen- 
tric distance of each GC system. Vrms.o and 7 are two parameters 
free to vary, with the latter representing the slope of the power-law: 
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increasing if 7 < and decreasing if 7 > 0. Uncertainties on this 
method were obtained by bootstrapping the sample 2000 times in 
order to obtain the 68 per cent confidence levels. 

Bearing in mind that the shape of the profile might change 
with radius (see the velocity dispersion profile of NGC 1407), in 
Figure [22] we show the best fit 7 for each GC subpopulation (in- 
cluding the combination of the two) as a function of the host galaxy 
magnitude. The overall slope of the Kms profile (7) of the the GC 
systems is generally positive, but it becomes < for brighter galax- 
ies, implying a constant or even increasing Kms profile. If we now 
consider the blue and red GC subpopulations separately, it is clear 
that this effect is caused only by the blue GC subpopulations in 
galaxies with Mk < —25 mag. The red GC subpopulations are 
always consistent with having 7 > and a linear fit to the data 
suggest no significant trend with the magnitude. It is not clear why 
we do not observe 7 < also for the red GC subpopulation. The 
two GC subpopulations might have different orbital configurations 
as a consequence of different formation histories and this would al- 
ter the shape of the Kms profiles. Therefore, it is important also to 
study the orbits of the GC systems. 

9.5 Kurtosis 

The galaxy anisotropy has classically been measured through 
the parameter j3 that gives the relative ratio of the tangential 
to the radial velocity dis persion for a given dynamical tracer 
jBinnev & T remaine 2008). The calculus of /3 for our GC systems 
is not a topic of this paper and we will instead keep things sim- 
ple by calculating the kurtosis k that measures the deviation of a 
velocity distribution from a Gaussian. In the regime of a constant 
velocity dispersion the kurtosis provides a rough estimate of the 
GC orbits. These will be will be radial (n > 0), isotropic (n = 0) 
or tangential (« < 0), respectively. Broadly speaking we measure 
the fourth-order velocity moment: 
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but in pra ctice we will use a mo re complicated bias-corrected for- 
mula from ljoanes & Gilll jl998l) . 

We divide the radial kurtosis profiles in two radial bins in order 
to have roughly the same number of GCs in each bin. In this case 
the required number of GC radial velocities per bin is larger be- 
cause the uncertainties are oc N 1 ^ 2 . We require at least 30 GCs per 
bin and this excludes smaller galaxies like NGC 1400 or NGC 821. 
As the calculation of k does not take the GC radial velocity uncer- 
tainties into account, a bias might be present due to the different 
velocity resolution of the literature and of our GC samples. Results 
are presented in Figure [23] 

We find that in the innermost bin most of the GC systems have 
kb < and kr < 0, suggesting isotropic or even tangential orbits. 
This is also true for kurtosis in the outer bin, but some galaxies 
(such as NGC 4365, NGC 4486 and NGC 4472) have the tendency 
of having more radially-biased orbits. In all cases the kurtosis of 
the blue and red GCs are consistent with isotropic orbits. 



10 DISCUSSION 

10.1 GC kinematic bimodality 

The origin of GC colour bimodality is still an area of debate. On 
the theoretical side, there is a lack of consensus regarding the lead- 



ing processes that drive the formation of GC systems (i.e. merg- 
ing, c ollapse, accretion) and whether one (see [Muratov & Gnedinl 
|2010|) or more of these processes (e.g jLee et alj|20ld) are needed 
to explain the observed colour bimodality. The contrasting numer- 
ical predictions are in part due to the diverse recipes with which 
simulations were built feekki et all2005Ll2008l:lMoore et alj2006l : 



iKravtsov & Gnedinl [20051 : iDiemand et al.l 120051)7 but they might 
also reflect a real diversity of formation mechanisms that built up 
the GC systems in the first place. 

On the observational side, the GC colour bimodality has been 
generally accepted because the two colour subpopulations differ in 
mean size and spatial d istribution, and they depend differen tly on 
host galaxy properties l lForbes et al.ll 19971 : IPeng e7aT]|2006h . Re- 
cently, the colour bimodality itself has been claimed to be an arti- 
fact of a str ongly non-linear colour-metallicity relation jYoon et all 
120061 1201 ]]), which would make the metallicity distribution of a GC 
system unimodal, skewed towards metal-poor metallicities. The 
resulting metallicity distribution function would be more metal- 
rich than the one classically obtained with linear, or broken linear , 
colour-metallicity relations jPeng et alj 1 20061 : iFaifer et al] 1 201 it 
ISinnott et alj201Cl : lAlves-Brito et alfcOllll 

We have found that GC systems around most of the surveyed 
galaxies have different kinematics between the blue and red sides of 
their GC colour distribution. This includes different kinematic posi- 
tion angle, rotation amplitudes and velocity dispersions. The kine- 
matics generally change smoothly with the GC colour, but we also 
observe cases in which there is a sharp kinematic transition corre- 
sponding to the blue-red dividing colour (see NGC 821, NGC 1407 
and NGC 3377 in Figure ll~2b. Our kinematic results are therefore 
strong evidence against colour bimodality simply being a conse- 
quence of a strong non-linear colour-metallicity relation. 

Most of the surveyed galaxies have small, yet significant, 
blue/red kinematic diversity (especially in their velocity dispersion 
profiles). The fact that we do not see any sharp kinematic transi- 
tion with colour in all galaxies is not surprising because these sys- 
tems lack intrinsic rotation of both GC subpopulations. Therefore, 
the rotation amplitude of GC systems can be insensitive for testing 
kinematic bimodality, especially for massive galaxies, for which 
one has instead to rely on other proxies such as the slope of the 
Kms profile (Figure[22l. 

10.2 The star-GC system connection 

Although the link between the red GCs and the host galaxy 
was somewhat expected from their similar spatial distributions, 
its kinematic confirmatio n was limited to a handful of galax- 



ies JSchuberth et al 



20 id : iNorris et alj|2012t iRomanowskv et all 



2009l : IStraderetalj201ll) . The results that we found for our galaxy 
sample are diverse in this regard. 

On one side we find excellent agreement between the kine- 
matics of the stars (and/or PNe) and that of the red GCs, indicating 
that the coupling between stars and GC systems holds also for less 
massive galaxies. On the other hand, we have also found galaxies 
in which the blue GCs behave as the galaxy stars or, more inter- 
estingly, as the PNe. NGC 3377 and NGC 821 are two examples. 
These are very similar field disk galaxies with a falling stellar rota- 
tion curve. However, the blue GCs in NGC 3377 do not rotate and 
are decoupled from the PNe. Conversely, in NGC 821 blue GCs 
mimic the minor axis rotation of the PNe, but not of the stars. What 
is most surprising in the latter case is indeed the fact that the PNe 
are decoup led from the s tars as well. Based on stellar population 
inferences. IProctor et all d2005h concluded that NGC 821 has un- 
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Figure 23. GC system velocity kurtosis. Symbols are as in Figure [T9l The plot compares the kurtosis k of the blue and red the GC systems in the inner (left 
panel) and outer radial bins (right panel). Dashed lines represent Gaussian distributions at k = 0. GC system orbits are tangential and radial if k < and 
k > respectively. 



dergone a minor merger or tidal interaction in the last 4 Gyr. Our 
results would support this scenario. However, it remains to be ex- 
plained why the PNe are not akin to the underlying stars and how 
this feature can influence the dark matter estimate of this galaxy 
dRomanowskv et al]l2003l) . 

How the star-GC system connection arose in the early times is 
still an open debate. Galax y mergers may con tribute to building the 
GC system of the remnant. iBekki et ai]d2002l) have shown that tidal 
shocks induced by galaxy merging can compress giant molecular 
clouds to form both new met al-rich GCs and stars . However, based 
on stellar population studies, iForbes et at] d2007t) have shown that 
most of the red GCs are consistent with being uniformly old, and 
few, if any, formed in late epoch gaseous mergers. 

A possible consensus scenario may arise from the idea that 
part of the red GCs might have formed dur ing the early turbu- 
lent phases of galaxy for mation at z ~ 2 l lKravtsov & Gnedinl 
120051 : IShapiro et al.ll201(i|) . At this redshift, galaxies are c harac- 
terised by fast rotating thick discs dElmegreen et al .1 [20091) frag- 
mented into super-star forming clumps. Observations have shown 
that th e origin of these clum ps is related to disc instability pro- 
cesses jWsn^s^^^]20Tll), but they may be also caused by cold 
flows dDekel et alj|2009h . It is thought that these gas-rich clumps 
migrate inwards to feed the nascent bulge (dragging along newly 
formed metal-rich GCs) and losing their rotational motion to dy- 
namical friction, finally acquiring a typical V TO t/cr < 1 for the 
bulge. During this migration, ~ 50% of the clump mass will be 
stripped off and it will end up forming a rotationally supported disk 
jBournaud et alj2007l) . 

This formation mechanism seems to qualitatively explain the 
spatial and kinematic similarity between stars and red GCs. How- 
ever, as far as we know, no numerical simulations have to date 
been produced to directly compare observations with the kine- 
matics of the red GCs forming during the 'gas-rich clump-driven 
phase'. Also, the st udy of dwarf ellipticals in the Virgo cluster 
dBeaslev et al 1 I2009T) and in the Sculptor group doisen et al ] l200l) 



suggested that disks may have been significant sites of GC forma- 
tion at early times. 

The poor radial overlap between studies of galaxy stars and 
GC systems complicates the kinematic com parison between thes e 
two. Extended stellar kinematic studies, a la lProctor et al. I d2009l) . 
will help to test the galaxy star-GC system connection. Also, it will 
be interesting to see what is the role of the PNe in this scenario and 
if the PNe-GCs connection in NGC 821 is the rule rather than the 
exception. 



10.3 Orbital anisotropy 

The GC orbital anisotropy, in this context represented by the kur- 
tosis ft, can provide other clues about the formation of GC bi- 
modality. Theoretical studies predict that dark matter haloes, in- 
cluding the baryonic trace rs within them, have radially-biased or- 
bits i n their outer region s dMcMillan etai]|2007l ; IPrieto & Gnedinl 
l2008t iDekel et alj|2005h . However, the small number of studies 
in this field, which have usually employed deep integrated stellar 
light and PNe, seems t o suggest that early- t ype galaxies harbour 



quasi-isotropic orbits (Gerhard et al 



de Lorenzi et alj|2009t iNapolitano et alj|2009l 1201 ll; iDeason et at] 



Douglas et at] |2007|: 



201 ll) . How does this compare to the GC anisotropy? 



Our analysis suggest that in the inner regions both GC sub- 
populations have isotropic or tangential orbits, whereas in the outer 
regions there is a hint that some red GC systems might have ra- 
dial orbits. We recognise that in both cases the GC orbits are still 
consistent with being isotropic. 

As already pointed out bv lRomanowskv et al] d2009h , the lack 
of radially-biased orbits for the blue GC subpopulation might be 
connected to the evolutionary history of the GC system within the 
galaxy. Processes such as evaporation, dynamical friction and the 
consequent disruption of GCs may change the anisotropy. The dis- 
ruption of GCs will be more efficient for those o n radial orbits be- 
cause they plunge deepest into the galaxy centre dBaumgardtl 19981 : 
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iBaumgardt & Makinoll2003h . Accordingly, they will become less 
luminous and eventually vanish altogether, making GCs on tangen- 
tial and isotropic orbits the most likely ones to survive and be ob- 
served today. 

The GC orbital anisotropy will need further in-depth investi- 
gation but the overall pattern suggests that, if the blue GCs formed 
from the accretion or minor merger of satellite galaxies, some dis- 
ruption effects have to be taken into account to explain the gen- 
eral absence of blue GC on radially-biased orbits. The scenario 
discussed above supports the idea that present day blue GCs were 
brought into the host galaxy via accretion and/or minor mergers. 

The key point here is that the different GC orbits mean that 
the resulting Vims profiles will be shaped by the orbits themselves 
(the mass-anisotropy degeneracy). Tangential anisotropy, for ex- 
ample, could mimic the presence of dark matter, whereas radial 
anisotropy could deplete it. Interestingly, we find that the blue GC 
subpopulations with tangentially-biased orbits tend to have increas- 
ing V rms profiles (suggesting a more dominant dark-matter halo), 
whereas the red GC subpopulation is consistent with having a de- 
creasing Vrms profile and more radial orbits. It will be interest- 
ing to test how a classic Jeans analysis that makes use of the dif- 
ferent GC anisotropics found in this paper compares with recent 
studies that do not require assumptions on the dark matter halo 
dWalker & Penarrubij|201 lh . 

10.4 NGC 7457: a case study for the formation of SOs 

GC systems have been exploited to study the unsolved conun- 
drum of the form ation of lenticular galaxies jBarr et alj [20071 ; 
lArnold et alj 1201 lh . In fact, it is still not clear whether SOs 
are quenched spirals: the remnant of a gentle gas shut off pro- 
cess that would preserve the kinematic signature of the progen- 
itor spiral <Bvrd & Valtonenlll990l : IWilliams et ai]|201Cl) , or very 
disky ellipticals: the prod ucts of violent merging events dBekkil 
1 19981 : ICretton et al.|[200lh that would diminish the original disk- 
rotation. The observational differences between these two scenar- 
ios are minimal and it also depends upon the galaxy environment 
dKormendv & Bendej|2012l) . SOs might still be rotationally sup- 
ported if they form in uneven- mass galaxy mergers, with mass-ratio 
ratio between 1:4.5 and 1:10 dBournaud et alj|2005l) . 

In this context, a particularly interesting case study is offered 
by NGC 7457 that we have discussed in $5.9\ This nearby lentic- 
ular galaxy contains in fact both signatures of a recent merger 
dChomiuk et al.1 [20081) and a disky GC kinemati cs with a possi- 
ble co unter-rotating core in the innermost regions dSil'chenko et al.1 
120021). In Figure l24l we compare the Kot/o" profile of the stars 
1 Simien & PrugnielluOOOf) and of our GC measurements (includ- 
ing the dataset of Chomiuk) with N-body simulations of merg- 
ing with different mas s-ratios of 1:3 dCretton et al .Il200lh and 1:10 
dBournaud et alj|2005t) . respectively. To do this, we select all the 
GCs within 1 arcmin from the photometric major axis and we con- 
struct two bins of 8 objects each. 

Results in Figuref24]show that our data are in good agreement 
with the 1:10 remnant, even though the discrepancies occurring at 
small and large radii seem to suggest that the mass-ratio involved in 
the formation of this galaxy was probably smaller than this value. 
These results, along with the counter-rotating core in the centre and 
the intermediate-age GCs, seem to favour the minor merger sce- 
nario. We also note that the SO galaxies simulated in B+05 are 
characterised by a strong rotation of the blue GC subpopulation if 
formed via merging, but we do not see any significant change of the 
GC kinematics with the GC colour in NGC 7457. This effect could 
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Figure 24. GC rotation rate in NGC 745 7. GC data from this 
work (supplemented with IChomiuk et alj |2008|) and long-slit data from 
ISimien & Pru gniel ( 2000) are shown as red and black points respectively. 
Also shown are the simulated merger remnant with a 1:10 (top dashed line) 
and 1:3 (bottom da shed line) mass ratio from [Bournaud et al] j2005h and 
ICretton et alj J200ll) . respectively. Typical uncertainties on the merger rem- 
nant profiles are 15% at 4 R e g. 

indeed be real, but it could also be driven by the low number statis- 
tics or by the large uncertainties (±0.2 mag) due to the conversion 
from HST to WIYN magnitudes. 



11 SUMMARY AND CONCLUSIONS 

We have examined the GC systems in twelve early-type galaxies 
(from lenticular to large ellipticals over a range of galaxy mass), 
nine of which published for the first time, with particular emphasis 
on the kinematics of their blue and red GC subpopulations. 

For the new presented data, we have used wide-field ground- 
based imaging (mainly from Subaru/Suprime-cam) and HST obser- 
vations, finding that eight out of nine galaxies have a significant bi- 
modal GC colour distribution. The study of the spatial distribution 
of these two subpopulations has revealed that red GCs are gener- 
ally more centrally concentrated and that they have a slope similar 
to the surface brightness of the underlying galaxy. This is in agree- 
ment with the idea that red and blue GC subpopulations trace the 
spheroidal component and the halo component of the host galaxies, 
respectively. 

Multi-object spectroscopy was performed for a bright sample 
of GCs using Keck/DEIMOS that provided an average velocity res- 
olution of ~ 15kms _1 (a factor of 3 better than most previous 
studies). We find a variety of kinematic profiles for both the blue 
and the red GC subpopulations. The salient results are: 

(i) The GC kinematics (rotation amplitude, velocity dispersion and 
rotation axis) varies with GC mean colour. In particular, we find 
GC systems which have a sharp kinematic transition at the blue-red 
dividing colour, implying that the GC colour bimodality is real and 
extends to kinematic bimodality. 

(ii) The rotation velocity and velocity dispersion of the red GC 



© 2012 RAS, MNRAS 000,[T]-?? 



30 Pota et al. 



subpopulation mimics the host galaxy stellar kinematics (including 
those of planetary nebulae). This property supports the scenario in 
which the red GCs form together with the bulk of stars of the host 
galaxy, for instance during the 'turbulent disk phase' at high red- 
shift. The rotation of the blue GCs is typically consistent with zero, 
but their velocity dispersion is always higher or equal to that of the 
red GCs, especially in the outer regions. 

(iii) The GC kinematics combined with other kinematic studies 
can reveal interesting features. We have found that the blue GCs 
in NGC 821 mimic the minor axis rotation of the PNe, but not that 
of the stars, suggesting that there might be a blue GC-PNe connec- 
tion that trace a recent merging event. Also, our data seem to be 
consistent with the idea that the SO galaxy NGC 7457 was formed 
via a minor merger with a 1:10 mass ratio. 

We have supplemented our dataset with ten additional GC 
systems from the literature and compared their kinematics to host 
galaxy properties. We have compared our results with numeri- 
cal simulations, finding no strong evidence that GC systems have 
formed via major-mergers. We find that the correlation between the 
Kms of the GC systems and the mass of the host galaxy (from 
lenticular to massive ellipticals) holds for less massive galaxies and 
that it is tighter for the red GC subpopulation. The blue GCs in 
more massive galaxies have increasing Kms profiles, whereas the 
red GCs have always shallower V tmB profiles. A study of the GC 
velocity kurtosis, suggests that blue GCs generally appear to be 
isotropic or tangentially-biased in the outer regions, unlike the red 
GCs which instead have more radially-biased orbits. 
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APPENDIX A: VELOCITY BIAS CORRECTION 

A feature of the kinematic modelling used in this work (equation[4]l 
is that the rotation amplitude is biased to higher values when the 
kinematic position angle is unconstrained. This effect becomes im- 
portant in systems with V TO t -C o, because the rotation (and its 
direction) is fully embedded in the dispersion. There are two ways 
to approach this issue: numerically, via a Monte Carlo simulation, 
or anal ytically, usi ng the Box's bias measure in non-linear model 
theory l lBoxlll97"ll) . The first method is the most straightforward, 
whereas the latter is comp utationally more com plicated. Here we 
adopt a similar approach to lStrader et al. I feOllI) and correct for the 
bias via a Monte Carlo simulation. 

In practice, for each data bin containing N GCs (see Figure 
lilt , we compute the best fit rotation by minimizing equation[3] We 
then generate 1000 artificial datasets, of the same size N, drawn 
from the best-fit model. We repeat the fit for each generated dataset 
using the original best-fit model as starting point. We define the 
velocity bias as the median difference between the computed and 
the simulated (known) rotation amplitude. This approach is adopted 
to correct for the bias also in Figure[T2] Tableland Table[7] 

We find that the magnitude of the bias correction depends both 
on the bin size and on the (Kot/o"). In the case of (V I0 t/u) > 
1, the bias is negligible regardless the bin size. For bins with 
(Kot/cr) ~ 0.4 (0.6) and N < 20, the bias can be up to ~ 
50 (40) km s" 1 . Within the same (Kot/c) values, the bias can de- 
crease down to ~ 3 0(20) km s -1 for N ~ 50. This agrees with the 
rule of the thumb of I Strader et alj d201 lh to identify the severity of 
the bias, i.e. the bias is small if (y vo t/v) > 0.55 x ^/20/N. If this 
condition is satisfied, we find that the bias is of the order of (or less 
then) our nominal median uncertainty of 15 km s~ . 
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